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1 1.0 INTRODUCTION 


Three computer programs were developed under the present work 
(see [1]). Two of these programs, PCASE and PCASLE, used the 
open wake models, whereas the third one, PCASLDW, used the 
double wake model of the partially cavitating cascade flow. 
PCASE requires the cavitation number, and PCASLE and PCASLDW 
require the cavity length as input data. The detailed explana- 
tions in regard to the type of model used for these programs 
are given in [1]. The general rule of thumb is that if the 
cavity length is short and/or the negative camber on the suc- 
tion side of the cascade is small, PCASE converges rapidly, 
whereas if the cavity length becomes long and/or the negative 
camber on the suction side of the cascade is large, PCASLE 
provides a better convergence. PCASLDW may be used when the 
pressure-side camber is large so that the calculated cavity 
boundary tends to intersect the upper surface of the blade in 
the supercavitating cascade flow configuration. : 


Six solution parameters are to be determined for PCASE: these 
include three transform coordinates, scaling factor for mapping, 
deflected flow angle at downstream infinity, and ratio of flow 
passage width at downstream infinity to that at upstream infinity. 
For PCASLE, one extra solution parameter, i.e., cavitation num- 
ber is added, making the total number of solution parameters 
seven. For PCASLDW, one of the solution parameters, i.e., 

ratio of flow passage at downstream infinity to upstream infin- 
ity, was deleted, but two new solution parameters, i.e., two 
transform coordinates, were introduced, making the total number 
of solution parameters eight. PCASLDW requires the most com- 
puter time of the three to carry out one iteration due to the 
increased total number of solution parameters. i 


In the following, the structure of the program including vari- 
ous subroutines, input data set-up, typical output data and 


listings of these three programs is described. Due to the 
similarity in the input data set-up of the programs, only 
that of PCASE is given completely. It is believed that there 
will be no difficulty in running PCASLE and PCASLDW once one 
becomes familiar with PCASE. 


EERE Ere eretryres 


2.0 STRUCTURE OF PCASE 


PCASE consists of a main program and several subroutines, brief 
descriptions of which will be given as follows: 


1) MAIN PROGRAM PCASE 


Specify the dimensions for data. 

Read input data. 

Exercise Newton's iterative procedure. 

Calculate lift and drag coefficients at the end of 
each iteration. 

o Calculate the cavity profile. 


0000 


2) SUBROUTINE OXFNEW(X,STOL,M,1I,DG,DF,FFF4) 


o Exercise Newton's iterative procedure in calculations 
for the six integral equations to find the six 
unknown solution parameters. 


x: Input and output data in array SXSI(I) 
SXSI(1): & - coordinate for the point B of the foil. 
SXSI(2): &€ - coordinate for the point C. 
SXSI(3): & - coordinate for the point F. 
SXSI(4): A, coefficient of the mapping function. 


SXSI(5) a,, downstream flow angle. 
SXSI(6): ¢€ = ratio of d, to dq, 
STOL: Control variables for the accuracy of Newton's iter- 
ations. 


M: Number of desired iterations for Newton's procedure. 
I: Counts the number of iterations of Newton's procedure. 
This is defined within OXFNEW. 


DG: This is one of the assigned finite differences for 
the numerical derivations of = However, it is 
no longer used in the calculations themselves as 
it has been replaced by the array DELI(I,J) which 


is read in at the beginning of the program. 
DF: Same as DG. 


FFF4: The residue of equation F(4); if FFF4 becomes larger 
than S4 of input data, the program is stopped. 


3) SUBROUTINE OFSIM1 (ANS,NOF,XCA) 


o Calculates integral I(1) of integral equation F(1) 
for special ease of foil shape with rounded left 
end. Called from subroutines: FILINTL, RMINT, CAVITY 
(see Reference [1] for F(1)). 


ANS: Final answer for the integral I(1) of equation F(1). 
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NOF: This is a controlling variable passed on from the 
calling subroutines: 


NOF = 0 - OFSIM1 called from FLINTL 

NOF = 1 - OFSIM1 called from RMINT for real part 

NOF = 2 - OFSIM1 called from RMINT for imaginary 
part 

NOF = 3 ~ OFSIM1 called from CAVITY OXFNEW at 
F(5) 


XCA: Integration variable passed on to OFSIM1 only if 
NOF = 3. 
4) SUBROUTINE OFSIM2 (ANS2) 


o Controls iterative procedure for calculating integral 
equation F(4). Called from OXFNEW (see Reference 
fad) for-F(4)). 


ANS2: Final answer of OFSIM2 
5) SUBROUTINE OFSIM3 (Y,XXII,IP,I) 


o Calculates g,(i) in integral equation F(4). Called 
from OFSIMZ. (see [1] for g,(é)). 

Y: Integration variable passed from OFSIM2, correspond- 
ing to é. 


XXII: Returns value of g, (&) to OFSIM2. The parameter 
is passed from oESIM2 to OFSIM3 in the form of 
one element of an array (XITC(I)) inside an 
iterative loop. 


IP: Number referring to the control point; IP = 1 to LPM. 


I: I = 2 for the subdivided middle point between the 
regular control points specified by IP; I = 3 for 
the control points. 


6) SUBROUTINE OFSIM5 (ANS5) 


o Calculates values of F(5) using Simpson's rule and 
Chebyshev-Gauss polynomials. 


ANS5: Value returned to loop in OXFNEW for equation F(5). 
7) SUBROUTINE FLINTL (YINT, KCTRL) 


o Calculates integrals in integral equation F(1l). 
Called from OXFNEW. 


YINT: Value returned for integral each time FLINTL is 
called. 


KCTRL: Control variable passed from OXFNEW directing which 
of the four integra).s in F(1) is to be calculated. 
(see [1]). 


{ 
1 
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8) 


9) 


10) 


11) 


SUBROUTINE G2 (XS2,AG2,IS2) 


o Calculates J» (€) in F(5) given integral variable 
€, i.e., XS2. Called from iterative loop in OFSIM5. 


XS2: Abscissa subdivision points from which g 
calculated, i.e., &. 


AG2: Value for Sy (2) returned to OFSIM5 after each 
time it is called. 


5 (5) are 


IS2: Number of control points on the second arc S$2. 


SUBROUTINE RMINT(SR,SM,MIQ) 
o Calculates rie toe rye 43 
tions F(2) and F(3) respectively. These values are 

used to calculate F(2) and F(3) in OXFNEW (see [1]). 


SR: Value for r returned to OXFNEW 
SM: Value for m returned to OXFNEW 


MIQ: Control variable passed from OXFNEW dictating which 
value (1, 2, 3 or 4) of r or mis to be calculated. 


SUBROUTINE CAVITY (XCC, YCC) 


o Calculates coordinates of points along cavity cross- 
section to give cavity shape. Passes cavity end- 
point coordinates back to OXFNEW. 


r mM, My, M3, M, of equa- 


XCC: Value returned to OXFNEW for x coordinate of cavity 
endpoint. 


YCC: Value returned to OXFNEW for y coordinate of cavity 
endpoint. 


SUBROUTINE IC2(SR,SM,XCA, ISIC) 


o When ISIC = 0 used to calculate ry and M4 of equa- 
tions F(2) and F(3) respectively. It is then 
called from OXFNEW. When ISIC = 1 it is used to 


calculate. 


SR,SM: When called from RMINT this is the returned value 


for ry and My+ When called from CAVITY, only SR 


is used and SM becomes dummy (see Reference [1]). 
XCA: Only used for ISIC = 1, integration variable. 


ISIC: This is a control variable which tells IC2 whether 
to do calculation for OXFNEW or for OFSIMS5 or 
CAVITY. 


0 called from RMINT. 
1 called from CAVITY IN OFSIMS for F(5). 
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12) SUBROUTINE MOSEC (A,B,ER1,X,J,XLPA,IS1I2) 


o Finds a root of f(x) = 0 where x must lie between 
A and B and f(A)>0, £(B) <0. 
A.B: A root of f(x) = 0 exists between A and B. 


ER1,ER2: Accuracy controlling variables where 


- x| ER1 and [f(x - f(x) | <ER2. 


Xyeal real) 


x: A root of f(x) = 0, found in this subroutine and 
returned to the calling program. 


J: Number of iterations done in MOSEC. 
13) FUNCTION AITKEN (XX,YY,X,N) 


o Interpolate the value corresponding to X with the 
data of XX(N), YY(N) specified by Aitken method. 


14) SUBROUTINE DETERM (A,N,D) 


o Calculates determinant of a matrix A of rank N 


A: Matrix input, requiring dimension. 
N: Rank of the matrix. 
D: Calculated determinant of A. 

15) SUBROUTINE ARCS2 (S2,XC,YC) 


o Calculates the arc length of the upper wetted 
portions $2. Called from OXFNEW in calculations 
for F(5) after the CAVITY subroutine. 


S2: returned arc length of are S2. 

XC: X-coordinate of cavity endpoint. 

YC: y-coordinate of cavity endpoint. 
16) SUBROUTINE ARCLEN (XSS,XL,XH,IS11I2) 


o Calculates arc length of small intervals between 
XL and XH along foil profile. 


XSS: Returned arc segment length. 
XL: Lower x coordinate of segment endpoint. 
XH: Upper x coordinate of segment endpoint. 


IS1I2: Control variable telling the routine whether the 
upper or lower edge of the foil is to be looked 
at; IS1I2 = 0 for the lower edge, IS1I2 = 1 for 
the upper edge. 


17) SUBROUTINE XCYC (XCB,YCB,CX,CY) 


o Calculates the point on the upper face of the foil 
corresponding to ine endpoint of the cavity. 


ae sane sa A 


XCB: X-coordinate of returned point on foil. 


YCB: Y-coordinate of returned point on foil. 
CX: X-coordinate of cavity endpoint. 
CY: Y-coordinate of cavity endpoint. 
18) SUBROUTINE BBBETA (XX,RBETA, IS1I2) 
o Calculates BETA(X(XSI)) 


XX: X-coordinate of the body for which the local body 
slope RBETA to be calculated. 


RBETA: Local body slope in radians calculated in this 
subroutine. 


IS1I2: Control variable; = 0 for the lower portion 
19) SUBROUTINE FARC (FAR, XLPA,X1B,IS1I12) 
o Calculates the difference between the arc length 
DSS and that corresponding the j-coordinates of 
XLPA and X1B. 
IS1I2: The same as that in BBBETA. 
20) SUBROUTINE SHAPE (X,Y,BETA,ISI1I2) 


o Calculates points along cross-section of foil to 
give shape of foil. Also gives the angle oi 
the tangent to the foil at each point. 


X: X-coordinate for which Y and BETA to be calculated. 
Y: Y-coordinate of calculated point. 
BETA: Angle of tangent to the foil at calculated point. 


IS1I2: Control variable to tell the subroutine whether to 
look at the upper or lower face of the foil. 


21) SUBROUTINE FC2 (T,F,XL,XH,IS1I2) 


o Calculates values of the function along the wetted 
arc to be integrated in the subroutine ARCLEN. 


T: Value to be calculated at. 
F: Value of the function. 
XL: Low limit of the integration. 
XH: Upper limit of the integration. 
IS1I2: Control variable; = 1 for sharp leading edge 


= 0 for round leading edge 


3.0 INPUT DATA 


The following data are those for the family program of PCAS 
(Partially Cavitating Cascade Cases) which include PCASE, 
PCASLE, and PCASLDW. Formatting examples are shown in Section 
3.1. Several data cards must be changed in order to run the 
different version of PCAS. 


It is important to note that these programs, particularly 
PCASLDW, were written to be able to handle blunt trailing 
edged foils properly. Readers are advised to see [3] for the 
cases having sharp trailing edges. The input data set-ups 
shown below are those for calculating the two-dimensional 
loadings for a partially cavitating propeller. The input 
parameters representing the propeller blade configurations 
include R, AAAA, to CCCC, A8 to D8 XROUND and A2AA to C2CC. 
The definition of these parameters is described in [1]. The 
thickness of the foil which was used for the plano~convex 
foil case in the previous project [1] is now a dummy input in 


these programs. 


9-14 


15 


16 


INPUT DATA FOR PCASE 


SYMBOL 


NGAUS 


TGAUS (T) 


WGAUS (T) 


XXM 


DELT (I,J) 


TH 


AAAA, BBBB,CCCC 


DESCRIPTION FORMAT 
Number of subdivisions used I10 
in Gaussian integration. 
Abscissas of Gaussian inte- 4F20.10 
gration. 
Weight factors of Gaussian 4F20.10 
integration. 
Weighting factor for solu- F10.8 


tion parameters in iterative 
procedure (0 to l). 


Increment for numerical cal- 6F10.8 
culations of partial deriva- 

tives. 

The thickness in percent of 3F20.10 


the plano-convex foil (dummy 
variable). 


End of the normalized foil 
= 1. 


Y coordinate of upper end of 
the normalized foil. 


Specifies the radial location 4F20.10 
on the propeller blade 

(normalized to be unity at 

the tip) 


Coefficients for terms in the 
equation of the cross- 
sectional shape of the lower 
face of the propeller blade. 
These coefficients are used 
in the second equation for 

x values along the cross- 
section where .2 < x < .8 
(see Reference 2 for the form 
of equation). 


a pecan aba ee Aa alan he nt 


re = 


17 A8,B8,C8,D8 


18 XROUND 


A2AA,B2BB,C2CC 


19 AAAAU, BBBBU,CCCCU 


20 A8U, B8U,C8U, D8U 


21 A2AAU, B2BBU,C2CCU 


22 IFLAG1 
NCHBY 
23 SBETA 
SBETA2 


Coefficients for third equa- 
tion of cross-sectional shape 
of the lower face of the blade 
where x < .8 (see [2] for the 
form of equation). 


Leading edge radius. This is 
actually used only when 
ISHARP = 1 (rounded leading 
edge). Otherwise it isa 
dummy variable. 


Coefficients for first equa- 
tion of cross-sectional shape 
of the lower face of the 
blade where x < .2 (see [2] 
for the form of equation). 


Coefficients for second equa- 
tion of cross-sectional shape 
of the upper face of the blade 
where .2 < x < .8 (see [2] for 
the form of equation). 


Coefficients for third equa- 
tion of cross-sectional shape 
of the upper face of blade 
where x > .8 (see [2] for 


the form of the equation). 


Coefficients for first equa- 
tion of cross-sectional shape 
of the upper face of blade 
where x < .2 (see [2] for 
the form of the equation). 


= 0 - for regular runs 

7 0 - for runs reading data 
from CASCLIM. Needs extra 

data for SXSI(2), SXSI(3). 


The number of Chebyshev-Gauss 
control points. 


Initial angle of incidence for 
a starting flat plate solu- 
tion in degrees. 


Body angle of a flat plate 


in degrees. Used as an 
initial solution. 


10 


4F20.10 


4F20.10 


4F20.10 


4F20.10 


4F20.10 


2I10 


5E14.7 


SF4 Always set = 10. Used to stop 
computation if the calculated 
are length Sl becomes larger 


than SF4. 
BETAB Body angle at point B. 
BETAC Body angle at point C. 618 


(initially assumed value) 


24 LPMS Number of control points 
over the § coordinates between 
& = -l] and b. Used for first 
a5 length Sl. (see Reference 
{1]). 


LPKS Number of subdivisions between 
= = b and the last point of 
! the coarse division made by 
LPMS. 


LPM2 Same as LPMS only used for 
calculations on second arc 
length S2. Note that there 
is only 1 segment spacing here. 


IFLAG = 1 - for first run which 
requires data to be fed in, 
i.e., but only SxSI(1) to | 
SXSI(5). i 


= 0 - for use of previous data i 
in which case data will be ; 
read either from a data card 

(1£ IREAD = 5) or from tape 

(if IREAD = 1). For IREAD = 

5, not only SXSI(1) to SXSI(5) 

but also SARC(I), BETAN(I): 

SARC2(I), BETA2 must be read 

from the data card. 


IREAD Used for controlling where data 
is read from. Either tape or 
card as above. 


ISHARP = 0 - for sharp leading edge. 
= 1 - for rounded leading edge. 


25 NITER Number of flow configurations 418 
to be calculated in 1 run. 


ll 


Number of iterations to 
stop the larger nest. 


Number of iterations for 


Newton's loop. 


Control index for varying 
either the set values of 
the angle of incidence, 
solidity, or cavitation 
number, depending on l, 
2, or 3, respectively for 
the NITER loop. 


26 ALFA1S Flow incidence angle in 4E14.7 
degrees (see Figure l). 


GAMMAS Cascade geometric stagger 
angle in degrees (see 
Figure 1) 


Solidity of the cascade 
( = °/s in Figure 1) 


Cavitation number 
= _ DL gpd 

= (By - By) /z0V) 

27 DE,DG,DF Finite differences for numer- 3E14.7 
ical derivations of = in sub- 


routine OXFNEW. These are 7 
replaced by DELT(I,J), no 
longer used. 


28 SXSI(I}, I = 1,6 This card is necessary only 6E13.6 
if IFLAG = 1; SXSI(I), 
I =_1,6 correspond to b, c, 
f, A, a2, and e. Values 
for SXSI(I) must be 
arbitrarily assumed and 
tried to see if a convergent 
solution is obtained. 


29 SARC(1), BETAN(1) Arc length vs. local inci- 2E14.7 
dence angle in radians for 
i the lower portion of the 
body; these data are 
29 + LPM needed only if IFLAG = 0 : 


SARC (LPM), BETAN (LPM) and IREAD = 5. 


30 + LPM 
SARC2 (1) BETAN2 (1) 


e 


30 + LPM + LPM2 
SARC2 (LPM2) BETAN2 (LPM2) 
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Arc length vs. local 
incidence angle in 
radians for the upper 
portion of the body; 
these data are needed 
only if IFPLAG = 0 and 
IREAD = 5. 


2E14.7 


3.2 INPUT DATA FOR PCASLE 


Several input cards change from those of PCASE. Only the 
changes are noted here. 


DATA 
CARD 
NO. SYMBOL DESCRIPTION FORMAT 
9-14 DELT (I,J) Increment for numerical 7F10.8 
calculations for partial 
derivatives. 
26 ALFA1S Flow incidence angle in 4E14.7 
degrees (see Figure l). 
GAMMAS Cascade geometric stagger 
angle in degrees (see 
Figure 1). 
SOI IS Solidity of the cascade 
( = ©/g in Figure 1). 
CAVLEN Specified cavity length 
( = 2, in Figure l). 
28 SXSI(I),I=1,7 This card is necessary only 7F10.7 


if IFLAG = 1; SXSI(I), I = 
1, 7 correspond to b, c, f, 
A, Ao" O, and E. "4 


’ 
| 


3.3 


INPUT DATA FOR PCASLDW 


Several input cards change from those of PCASE. 


changes are noted here. 


28 


SYMBOL 


DELT (I,J) 


ALFA1S 


GAMMAS 


SOLIS 


CAVLEN 


SXSI(I),I=1,8 


DESCRITPION 


Only the 


Increment for numerical calcu- 


lations for partial derivatives. 


Flow incidence angle in degrees 


(see Figure 2). 


Cascade geometric stagger angle 


in degrees (see Figure 2). 


Solidity of the cascade 
= C/s in Figure 2). 


Specified cavity length 


(= <5 


in Figure 2). 


This card is necessary only if 
IFLAG = 1; SXSI(I),I = 1,_ 


8 correspond to b, c, 


Oo, 


go, 


and h. 


15 


f, A, 


a 


2° 


FORMAT 


8F10.8 


4E14.7 


8F10.7 


3.4 TYPICAL DATA SET-UP FOR PCASE 


PIASCADSPTeTLC Ge CMI 200 es D4 STS MAT <hwa 
FETIAPS ePTASCK oL53eFRIPS 
FIT Is Se PLASC Ae AFTL eo FlIL 76 


e733 257222 02277835351 et7279E0S27 e3iCSE76 2 
oS$346 0825207 oe 7457335 6% oF RSLLESTIC eFLet46oFt 
eFetRTLIZ72 07931233332 
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4.0 OUTPUT DATA 


Most of the output are self-explanatory, however, those not 


explained in output data are described as follows: 


CAV (X) : 


CAV (Y): 


P(I,J): 


YINT4: 


SOLNR & SOLNM: 


F(X): 
SXSI(I): 
CLINF: 
CDINF : 
CCDD: 


CCLL: 


xCCcC: 


y¥Ccc: 


This is just a repetition of the input data 
TGAUS (I). 


Repetition of input data WGAUS(I). 
Solution parameters corresponding to SXSI(I). 


Each time these appear they are an updated 
version of those preceding them. 


This gives the x-coordinate of the cavity 
endpoint. 


This gives the y-coordinate of the cavity 
endpoint. 


Partial derivatives of Function F(I) used for 
Newton's method. 


Solution to 4th integral of equation F(1l). 


Intermediate calculated values of integrals, 
only used for checking the numerical accuracy. 


Residue of each function F(1)...F(8) 
Solutions. 


Lift coefficient at infinity. 
Drag coefficient at infinity. 
Drag coefficient. 
Lift coefficient. 
Cavity length to chord ratio. 


Arc length of the face of the foil. Either 
upper or lower face. 


x-coordinate of cavity endpoint. 


y-coordinate of cavity endpoint. 


4 


XS2D: 


XKSI: 


XXX2: 


SARC2: 


CP2: 


BBTAN2: 


XXX 2 


SARC: 


CP: 


BETAN: 


Intermediate values used as a check for 
progress of program. Can be ignored. 


Intermediate values used as a check for 
progress of program. Can be ignored. 


x-coordinate distance from the leading 
edge. 


Distance from the leading edge along the 
upper surface of foil at XXX2. 


Normalized pressure at the upper surface 
of foil at XXX2. 


Slope of the upper surface of the foil 
surface at XXK2. 


x-coordinate distance from the leading edge. 


Distance from the leading edge along the 
lower surface of foil at XXX. 


Normalized pressure at the lower surface 
of foil at XXX. 


Slope of the lower surface of the foil 
at XXX. 
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ARITESCS e955) AKL Le exT2exxtS eXXI Gel? 
AVITECSe5E) CXL Le KKT 2K KL Sg XX14 91? 
ARITS Chg SS) XML Le XXISo XXL So XXE Ge IP 
FOR) AWLwmmmae GC EL Ga Te 2X e2Ke 


lll 


i 


ange 


SIBIOITINE JFSIASCANSS) 
OLMENSLIV S2STCLOII Ee TAKCEACL IL) oc STC7)D 
COMMON YCOSeS3ETA2 
COMMON KITAC[IG) oe KI PVC 209d eS NSS5280299) oSA2526299) 
CIMMON CAVXC1LOO) ocAVVCLOOI oe FETASCSETAC a XCCCeNCAVeL PYMONS2 
SIAMIN AUCLIDMS TSCA eNCABV go ATTANCLION oFATANV2C19G) eFETANZE100) 
SOTHON FLAPANS IELTS gp JSAP OALF ALS GANS 
Ph hd | SSETA eK U4 oe ICPLy SAR6O9C513) 
CIATION [LOL eK K Zeki GF ANSe rete Vie IIL OS eK L ZL SS 3LGS eS 8 e0SS 
COMMON XSNCTI¢ TET OERCOVVY eKMeISERASSXSIOC 7) eSXSIOOCT) o YXS (7) 
SIMVIN PSTZeLAgSAICCSLS) oSAICICSL ID» OMe DE 
COMMON SETANCSISICALTAMCSIS2 eT Vel Me XITC200) ex JS(200) eXIX 
CIAGIN KRIIWDe A2ZAN 93233 0C252 
SIMMON BARNS BISSHCCIe A50539C39IBe TSAUSC100) oA GAIS C100) eNGAUS 
72a123019133926a4 
Tals SISIIITIVE CALLE FVII IEE de 
JSE SIMPSINTS RULE. 
33 1 IND = Le? 
1 eSTCIYO) = yxSCIMO) 
CIE = CISCIZTAD 
SJE = SINCJELTA) 
WS2 S42UL3 4avi A FASTIA OF Ge 
N3 25274945. 2742 
NS21 = NS2e1 
NS2a = VS2-1 
S2GAP = (XSTC32-XSTC222/NS2 
JS2 SF CISCALFALOSAN4 ASSIS CEST CSD HSA AISKXST OCT) 
33 2 [82 = leNS2t 
KS2 5 XSTC2S525A7%e(LS2-1) 
402 2 KS2*23¢ 
MYAS = XS2-XSTCHI#SIE 
KAYAS2 = XTAS@e? 
ASD = KSTCADeCOE 
ASO2 = ASJee?2 
PHIM SF OGAPEKAOSCOK FASS &ASO2) AKT) 
T= CIS2ceEdel)? 69 TO 3 
TF CLS2eE2eNS21) 33 £2 8 
CALL 32 (XS20ANSG201S2) 
32 CALCULATES 52 dIf4 XSI SIVENS 
232 5 IxXPCeANSS2) 
TiCT def 2040) ANSS2SC152) SaNSS2 
S2azRCIS2) = FS2eIdIK/II2 
39 TO 2. 
3 cavPiwe 
S2WERCL) = JWOX/SARECLetXSTCS)) 
ANSS2SCTS2=FALOSCSAITC Let KSE(S)I/ 552) 
33 T) 2 
4 CONTINJE 
S242 RCNS21) = Od9K/552 
ANSS2SCIS23 206 
2 SINPFINJE 
S2SRC1) = Oe 
33 10 JS2 = lLeNS2Ae2 
10 S2S%(JS202) = S$2S2¢552) 
LOCS2KERCIS2) HG eo @S2KTWCISSSLI MS2KCEVWISS[¢2eS25AFs30 
TFGlseNTe89) SI FI 93 
SARC2€1)=06 
33 50 LSAWE2eNS2e2 
50 S2SACISARC) FoF eCS2SICISARC1deS2SVUCC SARC Ol) 
32 30 ESA%lsi91S21 
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o 


ee een 


SC SARTACTSARII=SS2SRACIS ARS) 
90 CINTINJE 
ANSS = S2StCNS21) 
RITIRN 
Evo 
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SUBVOITINE LC2CSReS4ecSAeI SIC) 
OIMENSION «STC7) 
COMMON YClle SBETA2 : 
CIMMIN KITVC20G6) eK IF NC200) oANSS280299) opSARC26209) 
CIVMON TAVXCLIO ela VE LOO Me SZTASS3IET ACS ¥CIS a NTRS oi PMY ONSES 
SIAMOW AVCLODI eS ISAAI2ANCAS Vy 33 TANCLIO) »pBBTAV2C190) pSETAV2C100) 
SIMMIN FLMPANG IEPA ISM eA WAALS SANMS 
caso SBETAGKK4eICFHL e SAI5G9N6513) 
SOMMIN [DJn eg XAeKGoX Se TANG eT FeV oe VV eIHLSS oKLITTS IL GS oSMALS eDSS : 
COMMON XSNCTI¢6 SuLeERCOVVY oXMeITERAGSKSTOCT) opS¥SISIC7) eo VRS (7) 
SIM4IN FSLZ eh oSAILCS5L3) oSALSICEIL3I» PM eg dE 
SIMHOW FETANCS132 ee BE TAMC513) of Job AMX TIC 290) eK SIC 208) KIX 
SIMION KRISVI9A204932339eC25° 
SIAFIN ARRAS SIZZ gc rT AZeSZocA eGo FSASSCLIOI oI SEISCLIDIENGAVSE 
39 1 £9N = Le? c 
L xKSTCIONND & yXStTIN) 
MARL = XSTC&IYSINCEIELTA) 
YYt = CSTCO) #SISCIZF4A) 
YY1z vY¥lee? 
ti FIX MINT 
tw CAVITY JF JFSL45 FIV FCS AVEO IN TAIT 
STALLED FID4 FLINT FIR Feld.) 
FJ2 L2 JF F460 


Aa RO be 42 


aegaceue 


SR=06 
T4Zbe 
B4SxSTCLI 05 
2AMS=D3meXSTCz) 
347 Li3s4ele 
BS4MFs34-xSTC3) 
21155440734 
81223491/55 
313253445734 
TFQLSISeNE0e3)? 339 F3 23 
APLSECKC AML es WOXSTCLI on ch eC KCAMKSTOC3)) 
A722SK 2 Aen5STC2) 
A°SSS3RTCAF1/4F2) 

20 CINTINJE 
33.7 TSUMtlyNo43yY 
VASCAICESI4ISZ1L1) MCAT SIA) oLed 
VZSCAICTS IAI 43120908 COL5 54) 6313) 
SAS2SaRTCRASRI) 
SASTIAESAIT Lek L SIA) 9207893 
ASIPSIMN* AICI SIM) 3n 
X2XISKS TP KK] 
KPV¥IZEKOKIeeD 
WASKIAKOZevY12 
PIVSK2-CPA/SVI 2 
ASISYVL/SRV? 
T7OiSloeldel) AhReLelCKSTIVSA) 
TTCITSIC cE 202) RWR=a160 
Tee LST CeEde 3) AdeCLeHsSACZAFSD/( KSLF=xKSA) 
SR=2SKoSaseRaR 

7 SUSSN+SAZSeR aL 
2741230141532654 1 4 

! 


an 


SRESRePATSNCHBY 
SHSSMOIAT/NG ZY 
RETURN 

ND 
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iS a a a oth fitelnntah —«7 sate admndeerinietn* a 


SJSIQITINE FLINT YLT e<oTil) 

JE4MENSL IN K3STCO7) 

CIMMCNY YCCCeS3ETA2 

SIMMIN KRITMCZ[]9I) eK LE NVC 263) 9 ANSS28€290) pSARC26295) 

COMMON CAVXC LOO SCAVYOC 1590) s SETASe FETAC eXCCCeNCAV gL FMM VSZ 
CIMMIN ALSCLID) eLSAAT eo VOGT Ve IITAVCLIG) sSETANECLII) oSETANZ0100) 
COMMON FLAPANSDEL TAG ISAPOAL=AL OSAMA 

ca4toNy SBETAeXX4eTCIT 9 SARTOIC5S13) 

TISMIN LIZ eXM eX Se KT oTANGe rr ole VRoSFTGS oKL 31 5Se SFL GSoSMALSe OSS 
COMMON XSNC7)¢ TRE OERC OV VY a XMeITERAG SXSIGC7) pSXSIOOC7) o YXSL7) 
CIMFON PSEZeLMSARCCSLI) eSAWCICELZIoLAMDE 

COMMON SETANCS132 955 7AMCEL3) ole LO Ke XITC 209) oXSI0200) XIX 
CIMTIN KMRISVI9A244 9323395252 

COMMON AAMAS ASZSZI gT TTA eS eC3sISeFSAISC1OO) od SASSCLSOIeNGAUS 
SYSIGITINE FLINT CalLcULATis THE IVPEGRALS IN FC1) 

L548? SF 9 FIG SHAMV? Lele Il ese 

LSmARF = 1 FOR ROINIED Vole ILS 

TF FILLS AAVE BISVIED VeEes CHER VS HIV -GAUSS 

QJADRATUS 

SJAIRAT SRE FIAVILA TAN NOT SZ ISED*e SINCE SETA 

TS NOT A S4D0TS FINITL SNe 

NOASY = NIMSER OF CHEZSYSAEV-SAUSS ZIJADRATIZE SINTAQL FSINTS. 
PAL = 36191572554 

T2CTCLeé20e0) 59 TI 3 


eee oe 2 ce © | 


22°70 12 = 197 ° 
TO KSTCI2) = KSNOETQ) 
33 73 12 


9.99 11 IH = 197 
LL OXSTCI4) = YXSCIND 
L2 DONPIVIE 


Oo 


q S INL = OXSTC1 led t05 
IN2 = OXSECLIe1lede 05 
ALL = CON2*XSTC2)9/39N1 
ALS = CINZKSECSII/INI 

TS = OKSTCLIONS TQ) De 0S 

SMBS =ACXSTCQINKSTCL)) 905 
431 5 €3C5¢1L0e2/2435 
A3Z2 = CHBISeXSTCSII“IVSE 
FIAS = OXSTOZOKST C22) 205 
FO15 = CXSTC3+XSTO2)) #05 


AFL Ss CFCLS ele /FLAS 

AG2 5 CECLSSCSTOCLIIFF LAS 

SPACE2 = CXSTC3-KSTC20) 707 4M ‘ 
: VEAD 2744 FIR THE SITING AW 

BPCKCTRLeSEe2) SS TI 150 

LT? CIveSEe2) GO TI 139 

SCSPACS F Che@KXSTCLIIFFALIATC IPC) 

FSPALTE 2 TSPASE PEL IATL Od 8K) 

ta4si 

XCHOK = ele 

SPALTSISPATE 

39 20 ICHBY=1eNCHBY 

NOMEN AB YeITCH3 YO) 

AJCICAIVISTISC MSN H+o Le PAL C2eNCHBY)) 

MOST SIND SACI 543 4) #32 

PFCLTithWeldoi) 33 T) $33 

22 TECXC4CK este XKSI) GI TS 21 

TF CLO%seSEel PK) S2ACzZ = FS2ASE 

XCHOK 2 XIACK*SSATE 

134 = Tome. 
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33 TO 22 
a > oRSST EXISTS 37d KSICLO4"92) AND KSICIIM) 
21 TINFLINIE 
ISMA = [0%) 
Z3TANCICH3Y) = SETANCIIMDIOCZSETANCIIVIMAZETANCIIMAD? 
KoCx<STexCalxd/SPace 
' Py I3TAN TS JSED FIV CHSZvVCAdEv-SAISS INSTEAD OF ZETA. 
33 TO 20 
$39 Z3TANCICAHSY) = SBETA 
BITAN FIR TTEWWeEAel IS SPZTISELZI Lv IFSL%16 
20 CONTINUE 
120 SINTINJE 
TFCXCTRL ot ae4) 39 TI 4 
T= CMI TRLeta0e3) GI FI 3 
' IF (KITRLoTae2)? GO TD 2 
TF C1S4aR7elivei? 309 79 10 
YINS = Oe 
39 21219 ISJ" = 1sNCHSy 
A3ZD = CAUCISUMPFALLIZCAUCISI4) 412) 


110 VINT = YINT #@8TANCLSU4) ©SSRTCASS) 
YINT = YINTe#PAlL/Noa3¢ 
33 TO 109) 


‘ 16 SONTINJE 
os FACS LS FHF FASE JF AANVILINGS WVIED Le Fe 
NOF = 0 
XA = 
Chane IJFSIALCYINTe VO 9 XCAd 
re. KOA TS JUGT4Ye IJNoY JSZTI FIV F055) LNIKEVEde 
33 f3 1033 
2 SINTIWVE 
' RrA=060 
2 XSA 1S JUMMY. 
tsice 
CAwe LISS Ue SMeXCks ISLS) 
YINT=S2 
33 T) 109) 
3 CONTINUE 
woowot yTEStAL FIX de 
AJON) IS CALCULATES AND STORES 
YINT = Oo 
33.9129 ISJ4# = leNod3y7 
A331 = Lem4JISI4) 
A932 = COUCLSIVIGASLICCAZS2[©AICISI4VD?D 
SaA32 = SARTCAS2) 
A3C = A31/7S79A32 
129 YINT = YINT#aBCc 
YINT = VINE ePAI/ND437 
3) T2 1000 
eooooT VF IGtaAL FIR 14 
SUNVSE SETAC) STAN Tot AND rer age 
EXPECTED TI BE alLwAYS SMOIT4e USE SAUBS- 
SHES YSHEV SIJADWATIZE FIWJLhe 
AJON) IS ALREADY CALSILATEDs 
UF THIS 35 TRE FLAIST CASE ©9X ZLTAN26 
USE A SONSTANT ©92 JETAN2. 
B3TAN2 IS JSED FOR FS42vveSaISS INSTIAD OF SETANZe 
@ covfiuse 
DECI TERA esie2) 39 TI 150 
JF ClSeSie2) 62 FI 191 ¥ 
S6ETAZ MIST SE READ FOR THE FIRST WUNe 


eaarervandageeea gr 
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Rn 


33 183 Ic4+3v = Levcdayr 
130 S3TAVe Clroizy) = S3fra2 

NSZL=SNS2¢1 

39 195 LIsezLeVS21 
LE ZETANSCIOSISSBETAS 

32 T2 182 
150 CONTINUE 

Til Je3Ee2) GI TO 131 

To44 21 

xcecK = XSTC2) 

33.173 ICt3Y¥ = leNC43r 

, MASI = FCASeASOCI CARY) OF C15 

L352 LF OKC4lS4 eS 0exX<SI) 539 TI 151 

XCHSK = XCHCK * SPATE 

Li44 = 1044¢1 

32 T2 152 
131 CONTINUE 

£oe4a = L344e} 

BSTANQCIC45Y.) = SETAN2CIOM) 

POCBITANZACLIADHKSETAV ZCI IGWAD OOK SL eK SACKS SPAC 2 

ILMsICaByY 

XQST = FOASeAICIL Pe FSTS 

AVITTC5e250) Totes 3d Faneclitd eo x45 
256 FIRMATUL SN ge Toro l3Se oe SZSTAV 2S ee EL Ge Te 2K oeXXSlore i lees) 
173° CONTINSE 
131 CONTINUE 

YINT = Oe 

92 193 ISs4 = leNo43Y 

ADL = CSSTANACISUtPeFaT HCL eeASCISS4D) 

A392 = CASCLS IVI CASLI OM CAICISI49) 6442) 

$aaa2 = Sarvcas2) 
133 YEINT = VINE * 431753432 


LINE YINT#PAISNTASY 
L390 SINTINSE 

RETURN 

ZNO 
we 
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corneas Gas OER, 


SJSROITINE CAVITY (xSoeYCo) 
STAYS SUSIIITIVE TS CALLE2 FII4 IKMFNEw FI2 FCS). 
JIMENSION CKEXCLGI A SKEVCLIIPSANSIIC100) oS2I1206190) eS ICSI3¢100) 
DEMENSTIN SICAIGCL OID oe XSTC7) 
SLIMENSLIN CAVERCLIOD ecAv YY C100) 
COMMON YCOCeSBETA2 
TIMMIN MIP WCS90) KIT NC299)9 oe ANSS286290) ¢SARC206299) 
COMMON CAVXCLOODSCAVYC LOO SZETAS Se SETAC eXCOCeNTAd oL PUM ONS2 
CIWMDY AUCLODPsSESsA@reNVCASVeZSITANCLIIPS®SIBFAV206190) e3ETAN2CI00) 
SAMMI FLAPANG DELTA IZAP SALA AL GANA 
cua4oNn SBT TAoKXNMe TOOL es SAROO9C 513) 
SIMMDN LID oe KM OMSK FAN ser rat re VV eI SS eK STSS ec TLGS eS MALS gS OSS 
COMMON XSNC7)¢ ser eK RIO VVY go KMe I FERAGSKSIOC7) oSXSIOOC7 IO YXS(7) 
SIMISN FSLZeLAsSABOCSLIISSAWCICSLS I eLAMe DE 
COMMON SETANCS212) op 2ZTAMCS13) elie LAK eo XIIC 20S) oKJI0 200) go XOX 
COMMIN KRIGVIGAZAN 93253 —525°5 
SIMMON AARAG BIBS 9S ool y Ada s2e C350 dSeFFAISC1LOC) od SAIS CLOCI A NGAUS 
2 KTCS IS THE Cavity IND 2SINT CALCULATES IN SUSe CAVITY e 
wase * SISTIELTA) 
SOEL = SIVCDELTA) 
Pal = 30141592654 
02 1 .2AT1¢7 
2 eSTCLI4) F YXSCLQA) 
SSG4 = SGICLeecSTO5>) 
SSCL SAL OG C Lee KXST(SdI 40 2e#P Al) 
WoAd=36 
NCAVI=SNCAVOL 
SAVS 5 CXSTOC2POKSTCLII“SNEAY 
WZAVE THE JAST FIING QF “SI = 5 SINSE THERE IS A 
SIVNGILAGI TY FOR SLVSLE SAPLTUAL VITTEX 4IDlLe 
22 424 = IsCavi 
A = XSTC1) *#CAVS* (40Me1) 
AL FART QF JQAEGH = ZETAS WAT 
CKLMeEQe1> GO TO 3 
TFC<Mo Het aevVCAV1) 32 £2 10 
Towomnl(lS1CKSL) TALCILATION, COLLING JFSI4Le 
T= (LJeSE034) GI TI 75 


ie) 


gna 


oa 
maine? 


P] 
x 
z 
I 


NJF = 3 
SAce JSTHICANS eNOS ge XTAD 
= ANS 15 A SILIFION FIV TOLCKSL de KIT TS IGINTISAL FO xKohe 
TE (idveSGe27) ANSIIC<KL4Y) = ANS 
33 73 76 


75 avS = aNSI1¢K.M) 
73 SONTINJE 
Teens LI2CKSL) CALCSLATIONG 
T=] ye55034) 30 TD 77 


Ts1l2o 5 1 
CALL IC2CSReSMeXCAsIEIC) 
: ONLY SU TS VILLEZED|\= SH 15 FIL AIT. i 
TF CIs eEGe27) SRIZ6K.4) = SR 
32 79 78 


77 $® = STACK) 
78 CONTINUE 
sewoeetl3 OKST) TALSJLATI JVs S55 CHC3VSADVeGAUSS 
: QUADIRATURE FORMULA® 
3°CS = CXSTCLPeKST C29) 005 
S455 5 OXSTC22@XSTCLII eS 
asi C9Pl5%109/5%85 
432 GFF T5OKSTC3I0 05485 
<1 xCA*XSTC2) 


EM] = OXCM+L od eCKT ASKS TCL DOCK CMeKST C32) 
EK3 Ss SERTCEKI/SEK2) 

ZF33 2 O4950E<3 

TF Clse5e34) GI TI 39 

SICS = De 

3205 S94 = LeNoa3y 
EJLSCadCTSUM) A351) CAS 2-AUTT SUM) ) 

Sis1 = SOIPCESID 

TES = CleeAUCISUWI/SEI1 


SF3R Sf CHZSMAICISIYD O3F9C5OKTA 

§ SIDS = SITSCCEPSHSF3SzeSIUTCL eM AICISIMD eH 2d dF EF 5A 
SIC3S = SICS«PALISNENSY 
SISS 5 SLOTHALUSCOXSTC2P KOA EKCHMCSTELIIIOECS 
PiCluetae27) SIZSI3CK.4) = SITE 
32 fT) 81 

30 SIC? = SICSISCK.4) 

$1 CINTINJE 


secon ll aCKs ld owese 
SSE CHESBVSHiVeGAUSS JJASWATIZE FWA 
IN THE SAME MANNTT AS Paar FIR La CN 
QFSIM3. 


crpeogaee 


Li seSE034) 52 TI 32 
FOCS = CXSTCSP&XST( 23) #05 
FMC5 = OXSTC3SI$XST C20) e058 
AGL SF CFIT5 eT eI F455 

AG2 = CFPT5S<XSTCII FICS 
SIl4 = Ge 


33.7 ISUM=S LeNCHBY 
Ss CPETAVACISI4) +7 aL edrePAJCI SI) 
R322 CALCISIMD AGL HO CAICISI4) 6462) 
$35 = Surtcrp> 
= tassit3 
RO 2 FACSeAUCISUM) 42 255=xKr48 
= SISseRI72A0 
SIC4& = SIC4*#PAI/NEN3Y 
LiCl 2027) SITSIGC LI 4S SLS8 
33 To 33 
$2 SIC® = SISSI4¢< LM) 
33 SIVTINJE 
TE CLs eS Ge 27 eAND Kn Mer de2) ARITEC6055) ANSeSReSIC3sSIC4eA KLM 
TF CL set ee27 eA eho re 4O) AIF EC 5955) ANS eSToSTOSeSIS4 oKLY 
Te Clu etGo2 Feb eKh Me  GeS0) SRITECE, SED ANS eS%eSISZeSIC4 eKLM 
55 FIVMAL CLIK eg temelL lel 2el3el4 JF CAVITY AREssme me Gl Olde Te2XoeKxe 
Ae<udsee ITO) 
4CURSID 2 LEK ALREADY CALCIJLATEDe 
SS2 SF CISCKLFAL SOA SAAD TISCKST CSD OEAMMADSXST OTD 
SS FS CoANSSIALMOSROC SIS LAL OSC SJ2)794A Le SIC3 
LeSIS4/7 ALD /EK3 
329 TO 25 
3 325 2 J2TAZOP AL 
33 TO 25 
13 SSS3ETACeIAT . 
B2TA3 AVD JETASC 320Y AvS.ES AT 3 AND SP MUST BE S2EISIEED IY COMMONS 
25 CINTINJE 
KX4S = XCA@SIEL 
YYT = XCASKXSTC4O)#SOF. 
YYT2 = YyTee2 
Meu = CSTC9) 9 CDEL 
XXU2 5 XXJee? 
er3) = YyT2ecxJ2 


ue 


» 


oe 
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DdOK = OSASeKAS/ (XYZ OF AL) 


23¢ = 598¢3¢) 
S3C = SINGSS) 
S75 5 Jede/SCS4 


o4ieCKLM) = CSCHCES 
S<EY CAKLMD = SGCeCED 
2 fIvVlive 
CAVKXCLI=AD® 
SAVY YC LI =J6 
DI 18 IcAVv=SeNCAVie2 
TAUMKCI CAV) 5 CAVKXCLTAVHO2) CSAVS eC CCEXCLCAV$2) e408 
LOCKER GLCAVOLISSMEKCL SAIS 
1S CaVYYCICav) = CAVYYCISAV=2) 
LOTAVSe (SAD V CL SAVam2 Dore SKE CTS aval de SKEVCICAVIDS 30 
Li ¢IveESe27) SO TI 158 
33 T3 102 
133 39 102 ICAVZL—NoAvi4g2 
SAVMKRCEITAVISSAVXXCE SAV? 
L32 Savrqicavescavrr(icad) 
XCCOSIAVXONCAVID 
SSCstAVYCNSAV1) 
131 CONTINJE 
KTCSCAVKACNCAVI) 
YISSCAVVYONSAV1) 
ViTIK 
avd 


Wed 
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Sacrerge 


Cearesr 


‘a 


oneepecs 


SUBROUTINE S2 (xS208520TS2) 
JIMENSIIN KSTC7) 9 MT 2180239) oXI223 6200) oX125S0200) oX1248 (205) 


COMMON 
Satay 
SI44ON 
23449 
33443" 
COMMON 
~I440N 
c2449N 
ca44ON 
224408 
cosdoN 
fom Tk Fed] 


YCoCeSBETA2 

MIP90290) oC II NC290) eANSS280292) eSARC26203) 

SCAMCOCLIOPSCAMPCIOOISFELAZe SSS AS eK COC ey NIA piP V4 9 NS2 

AJOLOD) oT SAAT NCMIV 9 II TANCLIIO) sABTANZ6156) e3ETAN2Z€100) 

TL mAAN GIT oT Ae ISA OAL TALS SANSA 
SBETAGKA4M eg ICOL e SAWCOIC513) 

TOS eKMOXDeXT oe PANSe Peo Ve IZITSS eX F15S eo BL GS eS4A.Se0SS 

XSVC7)» CTO ERO GO VVV go XMe I FERAg SXSIOC7) ¢SKSIOOCT) eV¥XS(7) 

FST Ze uF eSABCCSL3) oSAICIC5135)9 9 P99 D= 

BETANCSLS oe SSCACCSLS) ol Joke ML LC 2002 96550200) OK IK 

MRIIVA9AZAL9032330C255 

BARR SIRS ero AD eT OCIA IGeFSAISCIOG) od SALISCLOOPe NGAUS 


TAIS SJSRIUTINE IS CALLED BY JFSINSe 
T4H15 SSSRIITINE CALSILATES FINSFISN S20KS2) atlin 


UvC.Je 


a32.25 


S$ L21¢6xXS2) T2 124CxS2)0 
XSI* AS2 IS [42 SILSFIIV JF CNTESRALS« 


32 1 La2F197 
L xSTCLardorxScrar) 


aalos 


301913592554 


COCLSAL OG Clee KSTIE)I SC 2eePAl? 
1 (IJe3f047) 35) 12 169 

e-eL 2101S) eeeee 
T4E SAME CVPESIATIIV AS T4ar IN 
SJSIISTINE TAVITY 72% 520430) 


VO = 3 
Tau. 3F 


xI21 2 


SIYLCANS eo VD % 9452) 
ANS 


12 (1Jyol9099) RIZISCIS2) = KL2t 
eee el22(0SL) an-20 

JS= THE SAE SUSRIUFIVE 122 AS 

JSE2 Av CAVETY dIT4 LSLCZL6 


Isics 


cal. It 


xlz2 = 


2037 eS 4 exS2eLSI5) 
SR 


NOTE TaaT S¥ LS JI*4y VARIAILZ® 
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SSIRIUTINE JETERY Cho Ned) 
> JETERY VEVISEG 02-26-75 
Van 4 

DIYENSTON 465095029 SAVEAC50050) 
LF OW efd0e 126) TD 95 


i 
309 9i st 
3 SAVEaCT 2) 
< 21 
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> RQe INTERSAANGE 
30 T3 23 
22 J esi ed 
23) SAVE<JS = SAVEACKeu) 
SAVERC Kos) = SAVZAC LJ) 
SAVEACL ed) = SAVEKJ 
£= €5 eNEe NWISO TI 22 
22 54 #} 
63 T 31 
30 fsztl+i 
M1 CONTINJE 
IF CSAVEACKe4) e&ae Bed $2 TI 28 C2-20-73 
4 = SAVEACT o<) / SAVZACKO 4) 
SAVESCI 94) = De 
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gu T2 35 
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RETURN 
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: ) TIMGIIN FLAPANG Ea PAe TSA 7 eA Ale SAMA 
' Cos4ON SBETA eKX4e ICAL SARZO00513) 
: CIMMON TOs a KA a KGa KS eg TANSo trea VR oISLSES oKLSTSSoSLGES eS4ALS DSS 
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XLPREXX 
33 73 21 
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COMMON 
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AJC1L0O? eI StAX 2 eg NCH3Y, S33TANC1I0) oSBTANZC100) eSETAN20100) 
FLAPANG DELTAS ISAP OAL ALO SAY4S 

SBETA gM M4 eo 1 S739S0R5096513) 
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PSU Zeb PF eSAVLCILS) pSHWCICSLS dei AMe DE 
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DIMENSION FLC20) pF IC269) CP QACLEL) oY XKSC 261) eFL2ACL O58) oe FGSCI 50) 
COMMON/FOILEND/XXOSeYYOS 
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COMMON /SCVTYL/CEVLENSEILSSE2 
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COMMON FLAPANeCEL TL ISAFeALFALOGAMMA 
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s2TTi€S—S551) SETAES 
eFiTElS e565) orhake 
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a TITECSeSE7) KMACUNG GA 7MhA ECE 
eITl Cre 23) AAbAL SS RESUNCrS 
eile Sl4) ACVeAcseleusuis 
eRITECSeE25) Acnade? 2325ueC25CU 

eRITECE go LS2E) LEEKS el MKS oSELT AG SPEAD a AN ChEY 

eM ITECSel%2Z4) Boot e I= oSF4 

et iTecSetS22) SSETAZ 

Fon aT Ola eg ebALALS eg FL reba c ka tec Seat eK Lie Sacv ar COCOUSHeF 1 ot) 
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° 
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FORM LTS o rg eR aeg He erg Xa tahhistorl  ehecXetsszstre= di eb gcxeo CCl Care 

ve 1 ceo? 

Ese FINYATC SIX eg MACS eo F lL ese c Xe rs Sree Pl Toby ck ae Ci ae gl ob aixetCcaeg lie 
Xs) 

S57 FURMATC So xg ae xROUN Ose gE Lebel Xe FAZAAS HEFL be CX eA SSE eer Lek ee ret 
KecCaeg = i205) 

TIE FOPMATCETIL)D 

AAAAGESEBSCCCCARE CONETANTS FOX 2eTIR4 CAMSERe ¥ AND SGPTCX) 

wwe m~wweeCAL CULATED FROVANOTHER PRIGREM CALLED * CAMS ER tw emmwwonn, 

499589C5 ANG DF ARE COEFFICIENTS FOR FOLYKRIMIALS FUP KX GPEATER TrAN eée 

£20923 AND CuD0K ARE Now JIMMY 

S*4 IS LSED FCR DETERMINING wHETHER TI CALCULATE SETAs 

1321 FOSMATCSEL4O67) 


O7098IO 


C IFLSG21 NEESS GATA CAROS FOR S¥SICl de Telefe IRTAD “MAY FE SET TS fe 
STF GF. AGet 6 SATA oILL SE READ TITHER FROM 

> DATA CARD TF TREaCSE 
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351 FORMAT CICIS) 
201 FORMATCGTIE) 
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So FIRMATCOE1L4 07) 
C CE eOGeOF ARE TRE INCREMENTS FIR CERIVATIVES I CKFNE we 
2 35tbeTo2 $ OFZLeE$S ARE JSEG SEFCREe 
225 SQRYATCIE1L4 67) 
1229 FOEMATCS Xe4hLFMel be 2X e4HL PK Re Loe 2Xeo GKSSETAZ ee lose 7 oEXe GIFT ACH eile 
R2kKqencmbYsegit) 
ZEEL FORMATCOOXe eGETAE AND BETAC AS FIRST GUESS=e 9 = li eS e2xeFlbo®) 
L324 FORMATCLOXe SHDEZeE 4 e 7 eo 2X oe SAD HOE 4s Te SHOFAGEl Go 79 fY oe GHSF Gz eE) 407) 
1521 FIRYMATCICKX es #SSTTA2Se9 E1407) 
SSETAZSSBHETAZ#PAI/1336 
BE TAB=SETABOPAL/I506 
BETACSBETAC#PAI/190~ 
5 oP AMSLOM2=\S2 
Le MMSLEM2 
NE 22242 
Lewes eewes 
eR ITECS 91489) LPM 20S 4ARF 
LGaF FLRMATCL UM e LPM Se gL 3g Qk gti StAX > He el 14 67) 
TSnARbst CDF CHAP Leke 
FIR FISNOTD wote 
SSETACPAT/1AL. 
TUKLZLehITIR 
VIZES FROM CrENE ae SUT IE THE LOSS SCES “CT GC THFCUTh 
4 OF P&ESET VALUE MUST 3€ JSiDe 
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Ta Ome ti me 
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ALFAIDSALFAIS 
SAMMASSGAMMAS 
SSLISSSCLIS 
TF CNh<eESel)? SO TI 24° 
LF CNMK eEGec) GO TS 241 
SOLIDSSCLIS+Gel*FLOATCIUKL=1) 
32 To 242 
241 GRMMADSSAMMAS*#2e#F LIATCLIUKL=1) 
G3 TO 243 
240 ALFALIDSALFALS<2e KF LUATCIUKL=1) 
242 CONTINUE 
x¥ayye 
AGFALSALFAID OP AT/1E5¢ 
Doh Sle/S2L10 
GAMMAZGAMMADOPAT/ Lele 
SELTASALFAL+GANMA 
FLAP ANZ Ce 
MRITE(CS 9655) ALFALDe SAMMAD SOU 
856 SSRYATILXGISHINCISENIE ANGLE=e 
XY¥39014%e07) 
antes e653) FLAPAN 
652 FORMATCSXeLIMFLAP ANGLE=eE1%07) 
STOLLZ20&=4 
STCLS=SeE=4 
ERC21eE=2 
CLES Le E 4 
wR ITEC eeS11) CAVLEN 
Bil FORMATCIGNe CAVITY CENGTHSeef1Ge7) 
C SPECIFY HYOROFOILTS CHAFACTERISTICS AND SEPe POINTSs 
XSsLe 
Y¥Csle 
X35=Ce 
XASle 
XMOOS20%SOCC 
YYODSASUSSSUSCEUSOBS 


C 
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CLG e TO LX GEMGAMMAZOE L407 9IX oSRSCLIUIT 
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ev erer ener op td erer ad 


as 


WRITECG e562) KA KGa AC oY Co KXO0s VYOO 


392 FIRNSATCI IX eSHCHORD Sel 14 e7e2X el TAUPPER SEPe PCINTHeE14e7e2xXe2CH COUN 
Ke PIINTOXCe YCOSC Elbe 7Toltoeciselel nds xXXODH* oF LC eb ecXe HVYC DS wy 


Y F106) 

START ITERATIVE PICCESURE. 
wwe~nGASIC FLOW IS THAT OF FLAT FLATEoremne 

TITERat IS INDEX FOR NUMBER OF ITER ATIINSe 
TTERA=1 
LE CIFLAGeE@e9) ITERAZ=2 
BIGS=06 
XAIGH=IO6 
XuOd S06 
TS112=0 
XINCRT=SXA/50 6 
00 248 LINC=14+50 
XwOw=axnIG4 
XA IGHESXLOW*XINCRT 
TALL ARCLENCSe XLOmwexXHIGHeIS1 12) 

249 SIGS=SI1GSe¢S 

eeeeFING 381GS2-<ee= 

FIRST CALL SHAPE TO FIND A CORRESPONDING TO CAVLENe 

XCCC=CAVLEN 
X*REESSCAVLEN 
CALL SHAPECCAVLENS YeSETAe3) 
wRITEC 592006) Y 

2030 FORMAT CSXe#VSrE1Ge7) 
YFREECZY 
yccc =yY 
CALL ARCS2¢(SIGS2eCAVLENG YD 
BRITE(S9504) 31GSe5135S2 

534 FORMAT (10K ¢SHBIGS=oE140795Xo e3TSS2S80F 1407) 

CALCULATION UF TRE BLAGE TRAILING EDGE THICKNESS 
CALL SHAPE C1 00 gUPENDS3SETAG3) 
CALL SHAPE (1 eOsONEND BETAS 9) 
THELSUPEND@ONEND 
WRITEC 502691) UPENOsSONENDs TIBL 


2031 FIRMATCIAXS®UPEND AT XSL meee te FS ober Xe *DNEND AT X=1 wmoweny 


X FSehe3Xe*BLADE THICKNESS AT XSL erwwnn gh 604) 

STOL=1eE=5 

wPMSLOMS 

LPKSLPKS 

LOM1=LP Ma} 

LPM3=LPM=3 : 
TlPL IS USEC FIR CONTROLLING FROGRAMS C FOR ITERe 1 FOR THE 
FIND XSITSeXSICeXSIF sAeAu=A2 BY JSIVS WVEMTONGS METIODe 
SxSICly=xSle 
SKSIC2)=xSIC 
SXSIC3)=xSI1F 
SKSTCG)SA weiCa IS THE COEFFTs JF MAPPING FCNe 
SKST CS =ALF2 
SXSI CSB 2=SIGMA 
SCSIC7TISXSIG CLOWER FIRST dAKE CLOSYER POINT) 
SKSICB=EXSIh CUPFER FIRST WAKE CLOSUER POINT) 

TF CTuKL o6£e2) GC TO 53¢ 

IF CIFLAGeEGe 0) GO TO 779 
TNATIAL GUESS FOR SXSI(I) IS cneeeoonwee 

READC5e768) CSXSICKK) oe MXS108) 
. GO TO 169 
779 READCIREADS G20) CSKSLCKX) @KxKH1 98) 
620 FGORMATC&F1I067) 
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$29 09 621 IC=leLPM 
621 READCIREAD| 622) SARCCIC) eBETANCIC) 
522 FORMAT(C2E1407) 
00 162i TCHLelPMM1 
1621 READCIREAD 2622) SARC2CIC) eBE TANZCIC) 
IFCIFLAG eb ae0) SO TI 480 
GO TG 481 
430 09 682 IST=1eLPM1 
432 BETAMCIBTI=HeSe(SETANCTSTISSETANCI3T#1)) 
431 CONTINUE 
150 IcprIsa 
sR ITECGs102) ITERA 
L322 FORMATO LICK eS LOAITERATION NOe=el2) 
00 65C inP=ieS& 
850 SXSIOCIRPI=ESXSICIRP) 
IF CITERAcGEe2) STOL=STOLS 
TFCITER ACES eMSTOP) STILASTOWL 


fomen ema 


CALL OXENERCSXST STIL oMAXITeIING OG ed Fo FFFS) 


aververner 


630 CONTINUE 

00 537 101=198 

MSNOTOLISSKSICIOI) 

537 wWRITE(6eS35) ICisSXSI(IC1) 
536 FORMATCIOX es SHSXSI( elle 2nd 2021407) 

CSPACE=Cle*SXSICLIIISEL SATO FAK) 

HCSPACZOeS *CSPACE 

FSPACE=CSPACE/FLOATCI2M@LOK) 

HF SPACHGeE#° SPACE 

ABETH—1L eo MCSPACE*FLCAT(LP Kl) 

cele. 
SPI=Q FOR FINGING SXSICI)D»g Tefee SXSICI9SYXSCII2 ICPIS1 FCP TRE RESTe 
ALCULATION GF PRESSURE JISTRIEJSTION PC>te 

IF CI TERA eZuGel) 306 TO 36 

DJ 25 lh=leLP4 
3& BETANOCISISBETANCIE) 

oJ 37 I5S16LPM1 
37 BETAMCCIBI=BETAM(I5) 

oC 385 IB=leLP MM) 
355 SETAG2tIa» = SETANZ¢18) 
36 CONTINUE 

UU 2=CIS CALFAL*GAMMAD STOS(SXS 105) 5A4MAD/ESPACE Cc 

D7CISXS/)AMMAGSIECISXSCSIC/AMMAGSIAFLACSOC=2UU 
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UCULESQRT C1 *SKS1(6)) 
@eULSUCULPCAVLEN*UU2* (1 eCAVLEN) 
Gal 2Z2QeUler2 
CP a= lew adl2 
IF CUTAUCED 1) CPwS-SXSI (CE) MCAVLEN OC 1 e=UU2Z2) #CLemCAVLEN) 
00 25 LG=1eL Pm 
wP=.G 

FIND CPOCASIP) WEXTe 

coeee FOR THE FIRST wETTED ARC PORTION Sleweee 
cP IS 3ASED ON UL AND Ple 
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CVAD PED OPEV HPEV CDE CV EDM area cr CV ever er dserererwg q000 


perverse 


LP=1 IS NEAR THE TeEe 
wIZLP4 IS NEAR THE Lefe 
IFC(LPeEQei) GO TD 521 
TF (LPeEGeL PM) GO TO 52 
S2SEK7 OKT TNCLAD) 


XITTNCI) IS CALCULATED IN OFSIM2 OF OXFNEW FOR FCO 6 


52 


521 
522 
25 


Q2=92"02 
CPCLP) SL e-UU22"Q2 
63 TS 522 

CP (LP) =-SxSI (6) 
60 TO 522 
CO(LPI=CP a 
CONTINUE 

CONTINUE 


EUNITNCC 

of=) MPL CFC 

YSeGALUK PCY oICXCNEKTIA=) LMP LCS 
DEeSLIX HCV eS CXCNEKTIA=) 22 MPL CFC 
ECAPSFedICISxS=321x 

ECAPSF to2@JICISXSSF91ix 

eC= 4CPCY 

YINOMPL CFCS D3CPCY 
YESMPLIFCH)V2CPCY 
DT4FLOPCHILCPCY 

YICISXS= 4 CP CX 

ECAPSFee2t) 2°0FCyv=3(PCX 
ECAPSFee2+d 1 CFC X= 20 PCX 

ELATSF we 7=V1CIEXS= DICPCK 

431 OT SG ILeQEARETICFI 
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eeeerneee tMAIN INSERT Leeeeeneaerrereneetnerererter ee akeee eter eteererareeewnent 


eeeeencP FIR THE SECOND ARC S2--e---- 
KUMBER OF CONTRC. POINTS ON S2 IS FIXED 
IN SUBRIUTINE JFSI4Se LeEos 
FOR BETA 


631 


€32 
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eerreevene MAN 
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1136 


HALF JF THE PIINT SED 
ANSGES IN COMMON = G2e 
92 686 NC> = LelL PM 
IFCNCPeEGel) GO TO €81 
T= ONS> eEGeL PMM1) GI TS 682 
G2 = EXPCANSS2S(NCP)) 
G2 2 a2*22 
C22CNIIFD = LemB2eUU22 
60 TO 669 
CR ACNCOP D=eSxXSI(é) 
30 TO 68C¢ 
CP2CNCP)=CPu 
CONTINUE 


AF OS ABS UFFFS) 

Te CaF SeGEeSF4) GO TO 1135 
33 TO 1134 

e@1TE(S91136) 


FOPRMATCEXseF (6) IS TOD LARGE TO CALCULATE BETA®) 
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GPeverarvereververer ar 
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FIND XXXCXSIP) FIRSTe 
21234 CONTINUE 
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erent 


T$1S2=9 


eceeFIRST SETA FOR ARC l-----0 


02 100 LLz=19L9% 

LPSLPM@LLPo} 

CALL SBBETACKXYXe8E T491S1S2) 

MXXXCLOVSHVN 

BETANCLP)=S5ETA . 

IFCLPeEGeL PM) SETABS3ZET 

TFCITERA LE SMSTOP1) SO TO 109 

BRITECSeAClL) uP oeSARC(KLF > XXX (LP) oCPKC LO) epBETANECL?) 
CONTINUE 

FORMATO L Xe 2nIoel Se lX eSHSARCH etl eT ol Xe MHXXXD gE LH eT elX eg Id4CPeEl daly 
XIX o6MBE TAN=9E1407) 


wew ee eMAIN INSERT Qet ee tae ae ee eect ee eK EEERE KEKE et RPT RH ETRE KHER EKER HD 
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meee eee eBETA FOR ARC S2eenen enn wen nn ne nnnnnene, 
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errs 


weeee 


CIN 


SARC2 HAS SEIN CALCULATED 
IN SJBROUTINE JFSIME ANDO 
STORED IN COMMON AREAe 


Tsis2 = 1 
02 429 LL>2=10.P4M1 
LP=LLP 


CALL 3B8BETACKXYX BE TAgIS1S2) 
TF (LP efGel)d GETACZBETA 
XXXZ2CLP) = XYX 
BETAN2CLP) = SETA 
TF CITZR ALE eMSTOP1) GO TO 329 
BVITE(Ge 229) LAP gSARC2ACLP) og XXX2(LP) pSOQCLP) pBETANQCLP) 
FORM ATCOXee Lowel Se lXe SARC 2=egbl4oe Zo 1X yg t@XXXZaee 
PEL Se Tel Xe tC P2=eoLl be Tol Xe tBETANQ=% oi l4e7) 
CONTINUE 
CONTINUE 
* 


wren a4 AInN INSERT Cervera eaeeaenreeererteeeeeetrertene tere arererenrenetennteet 


wewereWATN INSERT 3 tetera raree ee ee rere eeet rere ne ee eeeheeenee ener ke eetene 
* 
* 


D LIFT ANC ORAGe 


woceewewenF IT RST CL AND CJ FOR $1 PARTS 


1368 


USISO = SINCOELTA) 

usood = COSCJELTA) 

Ux8 = SxSI¢4)*UCOD 

UXE2 = UXBee2 

30 165 iTK = 19h PM 
LECT TK eSTeL PK) GO TO 106 

X>S = 1LeeCSPACE CFLIATCITK 91) 
6G TO 158 

MPS = XBSETOFSPACE®FLIATCITK=LPK) 
CONTINUE 

UXA = XPS@SXSI(4)#USTID 

JMAZ = UXhew? 
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eo “ oN aa Se “ere ileal 


PXXF = UCOD/CUXA2*UXB2) 
DdOX = OGAP OPXXOexPS/>aAl 
coeeTi = COSCBETANCIT<)) 
SISET1 = SINCBETANCITX)) 
OSIDX = SEXP CHXTINGIT<K 2) e Oe X/UI2Z2 
G1 IS CALCULATED AT JFSIM2 AS XITNCI)«© 
AND STORED IN COMMONe 
TECXPSelTede) OSIOX = -0S10x 
XLP1 = GS1DX*CPCITK) 
FLCITK) = @xXbLP1*COBET1 
SHCLTA = XLPLeSIBETL 
195 CONTINUE 
CeeeeeCL AND CD FOR S2 PART. 
NS21=NSZ+1 i 
NS2A=NS2-1 
SAP2 = CSKSTC32-SXSIC2) 2 /NS2 
32 °333 ITK = le NS21 
KAS2 = SXSIC2)¢GAP2eCITK=1) 
UXA = XRS2-SKXSI(4) # SID 


Lae oe) 


UKAQ = UXAee2 
PXXP = UCSO/CUXA2+UK32) 
CwOX = DGAP*PXXP#KRS2/PAI 


CSOEET2 ==COS(SETANZCIT<)) 
SIBET2 ="“SINCBETANZCITK)) 
DS20X = CxXPCHAKSGASCITK)) #Oe0X/5U22 
S 32 IS ALREADY CALCUAATED aT GFSIMS AS 
S$ ANSG2S(I)e STORES IN COMMON AREA 
AXLP2 = OS20x#CP2¢ITK) 
FL2CITK) sexLP2*COBET2 
FIQCITK) = xXEPeeSI3=T2 
336 CONTINUE 
SPACE = CSPACE 
CAIFT = OeS*#CSPACE FFL (2) 40 eS SPALTE®FL(L PMI) 
CORAG = Ged*CSPACE 8F3(2)4005*© SPACE RFD (LPM1) 
DSO 11 IVA = 2etPM3e2 
TFC TUaeGEeL PK) SPACE = FSPATE 
CaIFT = COIFT*SFACE SC FLO IUA) Ho etl CIUASIIGFL OT UAt2II/ 36 
211 CORAS = CORTAGHSPACE*CFOCIUA) t4erFI (CI JAFL FFI CIUAS2)I/3 6 
53 321 IVA = LeNS2Ae2 
CEIFT = CALF T#GAP2 eC FL 2CITUAD +4 ot FL AC IUATIIOFLZCIUAS2) /26 
321 CORAG = CORAG+GAP2Ze(FD2CITUAD H4e#FO2CIUAS1IMFDSCIUASZ) /te 
erereraDG THE FORCES ON CAVITY PORTIINSs 
SUPROUTINE xCYC CALCULATES 
THE POINT OV THE UPPER SLADE PORTION SORRESPe TI THE CAVITY ENC POINT. 
CX ASXCCC 
Cryazyvccc 
CALL KCYCOXCCCBe YCCCS eC xXAgCYA) 
CLIFT = CLIFT*#SxSI(&) exCCCB 
CORAG = CDORAG@SXS1(5)*YCCCB 
wee eeASUUST CORAG FOR THE 3ASE PRESSURE We 
CORAGSCDRAG@CPw* THEY 
\eceeexCOC ANG YCCC ARE T4E END 2PDINTS OF CAVITY» CALCULATED IN 
SUBROUTINE CAVITY 
STOVES IN COMMONS 


42690 


Oo 


e 
* 


eeeeeeereeMA IN INSERT 2 tee eee eee eee eet Se eEREHHEHAHAEEEER HEHE TERE AERA EEREEEES 


Cvarvaererar,m cranes 


FIND BINF IN 2e1le 
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UZULECOSCALFAL GAMMA) /SCOSCSK¥SI (5S) GAMMA) /ESPACE 
OC aN=COSCALFAL*GAMMA) #COS<(SXSI(S)*+GAMMA) 
BINFSOeS#SINCALFALSS XSI (65) %202GAMMA) /DOUN 
BINFOATAN(CLe/SINF) 
AINFSOe5*PAI@BINFOGAMMA 
oc COSTAR AND ALSTAR ARE BASED ON VELOCITY AT UPSTREAM INFINITY IN (x oY) 6 
COSTAR=CORAG 
CLSTAR=CLIFT 
UINFSO0eS#SQRTCLetU2U1L #8242 e#U2U1 #COS CALF AL=SXSI(5) 9) 
FINF S20 eO0GAPeSINCALFALMSKXSI CS) /CUINF *COS(SXSI (5) *GAMMA) ) 
CLINF=CLSTAR*CCSCAINF) =COSTAR®SINC AINE) 
COINFSCLSTAR*SINCAINE) *CDOSTAR® COSC AINF? 
CLINFSCLINE/UINF #22 
COINFSCOINF/UINF ee 2 
ARITECS e117) CLINFeCOINE 
117 FORMATCIXe34¢HCLINE OR COINFSFORCE/I/ 22COUINF ## 29 5X eOHCLINFHeEl 4e7 
XIX eSHCDINS 921407) 
SRITEC5 9113) FINE 

LLE FORMATCIX s34HFINF IS OBTAINED FROM “MOMENTUM EGNe EHFINF=eE 1467) 
wRITE(6 9221) 

221 FORMATCIXe4EH=-"CCLL $ CCDD ARE BASED ON U1 IN ALFA! DIRE e==-) 
COLL ECLSTAR*COSCALFALIMCOSTARESINCALFAL)D 
CODOSCLSTAR*SINCALFA1) *COSTAR®COSCAL FAL) 

&LOGSCCLL/CCOD 
PR ITECG L321) CCODe Corr eALODeU2uU1 
152 FOPMATCLXeSHCCOD=e E14 eT elLXeSHCCLLig= Lhe 7 ol Xe4nlL/D=H ek 14 eT ei Xe 
Xx SAUQU1=0E14e7) 
MSTOPISMSTOP=1 
TFCLTERASLEOCMSTOP1) 32 TO 140 


seetewter a4 Aly INSERT 6 terme a eat teeter ee ERE EERE ER EERE REREKEREREUE EPEAT ERE ee 
* 
* 
CAVITY SHAPE 
ALREADY CALCULATED IN 
SUBRCUTIVE CAVITYe 
WRITE(6 9257) 
237 FIRMATC Eke reaCAViTY SHAPE omewnn) 


vos nasanneh te Sin Ae nee At es malt ee A es ee 


moa danoan 


NCAVLSNCAVL1 
3S 265 KCAV=1eNCAVig2 
23 WRITE(Ee 236) CAVXCKCAV )eCAVYCKCAV D { 
235 FIRMATLLIXe *XEesEL4e7el OX #VEeel 1467) j 
: * 
; Tewewerenead Aly INSERT &© seeeeretemeenerer ee errektee ree Reee ste eeeet eee ewer ee 
140 CCNTINUE 
XCCESD6 | 
YOCC=96 3 


wR ITE(6 0823) 
B22 FORMATO sg temeeewe meee J POER BIBY SHAPE ce eemnmnne) 
3G 621 ISHP21¢e81 
X3 002" (iS4P-1) 
CALL SHAPE (Xe YeBETAs2) 
921 wRITECG 9522) XoY 
D22 FORMATUS Xe *XSeoF 10 eS e2X oe VeroFl 005) 
REaINO 7 
WRITEC7T 9763) CSXSI CK) @KXE1 08) 
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756 
737 


17S6 


33 


vw 
ow 


FORMATCSF1Ce7) 
09 766 IC=1lslLPM 

wRITEC70767) SARCCIC) oBETANCIC) 
FIRMATCZE14 07) 

05 1765 IC=leLPMM) 

aRITEC7Ts 767) SARC2ZCIC? eBETAN2CIC) 
LFCITER AGE eMSTOP) 33 TO 3999 

LPKISLPK-1 

SPACESCSPACE 

HSPACESRCSPAC 

00 50 IM=leLPM1 

TF CIMeE2el)? GO TO Si 

TIF CIMeEGeLPM1) GO TO 55 

TF CIMeEGeL PK1) GO T3 97 

TF CIMeEZGeLPK) GO TOD 99 

TFCIMeGTeL 2K) GI TO 92 
KXY¥S—LeeSPACE#FLCATC IML +HSPACE 
XZCLIS—Le*SPACE*FLOATCIM<2) 
XZ2(292K2012¢SPACE 

XZC33SKZ02)4SP ACE 

XZ(4)=KX2032 SPACE 

60 TO 99 

SPACE=F SPACE 

nSPACESHFSPAC 

XYSKSET#MSPACE® SPACE*FLOAT(C IM=LFK) 

RZ CLISXSETOSPACE wFLIATCIM@L2K=1) 
M2Z2€2)2XZ01)4SPACE 

XZ032=KZ(2)*SPACE 

XZ204)2X203)4SPACE 

00 56 IK=194 

YBECIKISBETANCIM*I K=2) 
BETAMCIMISAITKENCX Ze YEE OXY 93? 

329 TO 161 

BETAMCLPK1) =005*(BETANCLPX1) 49 ETANCLPK)) 
GO TO 181 

BE TAMCLPKI SC eS #(BSETANCLEK) *SETANC(LPX 41)? 
33 TO 151 

BETAMCL S005 *(BETANCL) *BETANC2)) 

669 TO 151 

SETAMCLEM1 PSOeS#(SETANCLFEM1) *3ETANC LPM)? 
CONTINUE 

CONTINUE 

TFCITERASEGe1) GC TO 6 

DD 41 TE=1e.PM 
BETANCTEISBETANCIED # CL e=XXM) *BETANOCIED @XXM 
060 42 IFG=1l_5LPM1 
SETAMCIFGISSETAMCIFS) #01 ee XXMDSBETAMOCIFG) xx 
09 425 IFG=leLPMM1 

BETANZCIFG) = BETANQZCIFG) #C LemKXX4XDOBETACZCIFS 2 #XXV 
02 852 IRP=1e8 
SxSICIRPI=ASXSI CIR # C1 ee KKM SSXSTOCIRP De XXM 
ITER ASITEVASL 

IFCITERAeGTeMSTOP) GO TO 286 

3) TO 160 

eRITE(6929) 

FORMATCSXe26HITERATIIN WAS TERMAINATED o) 
CONTINUE 

STOP 

ENO 


SJUBROUTINE IKXFVERCXeSTIL ge Mole I5eD% 07 FFG) 
DIMENSION £08) 9205003) 0X08) 9 QC 508) oe XRRICG) Oo XMITES)D 
SIYMON/SELTAGSOZLT(3 93) 
COMMON /SCUTYLICAVIE VS 31SS2 
COMMON/FREECAV/XFREE Ce YFREECS 
ZIMAON YCASeSZETA2 
COMMON XIT#C 200) oXITNC 200) eANSG2SC200) oSARC2 (2003 
COMMON CAVXCLO0) CAV YC1O00) eZETAS SETAC oXCCCoNI Ad op FUNG NS2 
COMMON AVC100) pI SHAR» NCHBY, 8B TANC1I20) pBBTANZC100) eBETANZ C100) 
CIMMYT FLAPANGIJELTAeISAP eA FAL eS GAMMA 
CIM4oN SBLTAgK C49 ICI, SARC096513) 
COMMON LOULe XAaXGaXleTANGelP oCoYReIBIGS oXLEISSeSIGSySMALSeCOSS 
COMMON KXSNC3)D9 SLOT Uo VV eXMeLTERAGSXSIOC 8) oSKXSTISCE) oVXS (8) 
COMMON PSLZeLP pSARCO513) eSAWCICSISI¢ LPMe DE 
COMMON GETANCS13) eBZFAMCS1L3) ol el Kg XIT C260) eK II 6200) 0 XIX 
COMMON XRIUND yg AZAR 93238 eC2C0 
ZIMGON AABN DZ DB eC CDT AS a TIeCIeIGe TSAUSCLOO) eA SGAUSCLONI NGAUS 
CINMHINV/SPC/ESIATE 
COMMON/ TAUL/ TAU 
Pal= 30141532652 
Ts0 
TPCITEUA ce. E03) GI TI 272 
33°57 Llv=les 
57 atITE66e65) Ll ue XCILU) 
55 FIRVAPCIX es 24K eT le 2d sel lde7) 
272 CONTINUE 
Sil 7seLer2et JETT 791) 
3& Slis2e*DE 
3l6220906 
IF (xO? eSTeSI7) XC7ISSI7 
LECxCL eb TeSi1) KCLI=ESTA 
SILOS" C1) 420896 
TOxC2) ok TeSI10)2 X¢202S110 
SILL SExnC2.¢2e%06 
T& (x03) oh TeoST11) *C322ST11 
SI82xC32*2ee OIL T(3el) 
T= (x3) ec. Te SI3) XC322S13 
TE GxCSd ol TeSIS) xX(4d2ST5 
SISSC0 eS e2AT GAMMA # C1 0002) 
ZPOCXCS dol Te0e) GO TO 7h 
TFOXC5 eGTeSI5) xCS)S5T5 
GO TO 79 
Te TF CasSCxC3  PeGTeSI5) nCSdSeS15 
73 CONTINUE 
T2QxCod eb feeO01) X0525e001 
39 56 [1 J=108 
$6 WRITE(6966) Il Ue xCIIu) 
Tys1 
Sewwaanlt ( L) cre see cme ewes ce cceseewe econ wccs cm menerraee woes noes ne em owwe enn 
33 20 IK=198 
20 YKXSCI4) 2x14) 
§ CUNTIVUE 
€oTRL SF 1 
CALL FLINTLOYINTIsKCTRL) 
: SJBVOITING FLINTL CALCULATES THE INTEGRALS IN FCld6 
42TRL = 2 
Sale FLINT. CYINT2e425TRL) 
Mort = 3 
TALL FLINTO CYINTS e<STRLD 
“CTU Ss 4 


70 


320 


lo tarantino 


CALL FLINTL CYINT&sXCTRL) 
<oTRL=5 
CALL SLINTUCYINTS eKCOTRLD 
<cTauss 
CALL FLINTLCVINTGe<CFRL) 
COCLSALOG CL o+V¥XSC6))/020*PAl) 
CSIFALIGKCCISCYKS CS IHSAMMAISSISCALZAAL ©GAMMA DBESPACE) 
SSG=SURTC1Let+YXSC6)) 
J2SCISCALFALSGAMMAI/CISCYXS(5) e¢GAN44) SESPACE 
U22=U2e22 
J22TSL0/U22 
THOSSSS*CAVLEN/UZ¢ (1 eo CAVLEN) 
TAUs DZALOS( TWO) 
Te bsSAT Lem Le~VXSC5)°UZ2I) #CAVLEND 
TAUBI=ALOGCT HL) 
TF (UTA et S00) TAU WsTAUd 
TF CUTAUseE Gel) TAUWSTAUML 
FA S SCYVINTISPAIOYVINT2°¢CCOCLOSSI/PALD *YINTS 
LOVINTS/PALWYXSUS TASS MCVINGIS@VYINTS)/27Al) 
T= (IveEGel) WRITE 65079) YINTLeVINT2e VINT3eVINTS 
FIRMAT C1OKewwemmlTi gl 2ZelZel4 GF =CL) AREwmm ate HC EL Ge Te2XK)) 
TE ClueEGel) FC1) = FA 
T* (1Je£Ge2) 60 TS 3 
IF (Ide Qe3? SO TD 4% 
= ClJeEGe4)? 50 TO 329 
F (1veEG05) GO TD 321 
F (IlJeEGe6) 60 TO 322 
F(TUebSe66) 69 TI 3222 
IF CIUeEQe460) GO TD 400 
TF (Tl JseEGe401) GO TD 4021 
F CL sebae802) GO TI 402 
TE CldeEae 443) SO TO 443 
TYS=TAVCYXS(S) *¢5A%4A) 
DT SDAS“TYE®CAVLEN® SSS/THO/U2 
TS CuTAveZ del) JTIASOCAVIENETYS/J22/¢0 T1982) 
OTBISzeS/IZeCAVLEN/SSS/TWO 
T= CUTAUsES el)? ITHOSHPeS#CAVLEN/ (Td *%2) 
POL95) = TAVNCYXSCSI CSAMMAD A VINTS SIAL or ee CVINTS Mf INTSI/PAIOCTRDA 
PC Le OIS—VINT 3/C2etPOl eC LetV¥XSCE) DI FEYINTS=YINTS) SPAT OO TWDS 
TJ = 2 
YKXSC1) = XCLI*ODELTC1¢1) 
32:79 5 
FIP = FA 
2 ) 
YXSC1) S KCLISDOLT (lel? 
Go TO 5 


PQ1le1) = CFELPOFIGI/C2e*DELTC1:e1)) 


x1) 
KXC22*DELTOC192) 


FIP = @FA 

YxSC2) = KC20-OELT C102) 
bd = 5 

32 f3 5 

Fi@ = «FA 

Ptls2) 
YxuS¢2) 
YxS¢3) 


CELPOFIG) SC 2e* DEL TC 292)? 
KC2) 
XC3ISOELTC193) 


oto 
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Iv 26 
69 TO 5 
322 =19 = -Fa 
TU=$6 
YXSCSVSKCSI-DELT(193) 
32 T2 5 
3222 FlQ=-FaA 
YXS03)=xC€3) 
POl93) = CFIPOFIOI/C2eeCELT(193)) 
Ple4) = De 
TJ=403 
ADSX7SABSCXC TID 
YKSO7ISXC7TISCDILT CL eZ #ASSKT 
60 TO & 
430 FIPs-FAa 
1J=401 
YXSOTISNXCTIMPELT (197) ©a3SX7 
$2 fra 5 
431 FIGS-FA 
PCLe 7ISCELIWMFIQI SC SoM DELT C1 e7) tABSX7) 
YXSO7)=X€7) 
TJ=4 02 
i YXSC8=NC3) ODELTC1 03) 
; 332 72 5 
492 FL>Is~FA 
TJ2443 
YKXSOCSISXC3IV@DELT C193) 
69 TO 5 
a3 FLQs=FA 
PCLSSISCFLISFLAIM/C2eeOEL T4158)? 
YXS¢€8)=xX(3) 


Ww 27 
380 CINTINJE 
XKKXZALOG(COSCALFAL*SAMMAI/TOSCYXS C5 )+GAMYA)/ES2ACE) 


XKX1 2 YXSCG)*#SINCOE.TAD 

YYL = YXSCOd#COSCIE.TAD 

Yr lzsvvYiee2 
COCLISALOGC Let YXSC6II/C2eeP AT) 
ToL = COTLMSXKAKX/FAT 

XRR = Oe 

X44 = Oe 

33.331 “ID = lesb 

CAL. UMINTCSOLNRe SOL V4 se MIG) 
MIRTCALS) = SILVr 

XEMI ¢41Q) = SOLNM 

MIRI Ss HKAVICAMIQI/ PAL 

MMH = mKXMYTOMIQD/24T 

LF (ML AeEDeL) KRARZ 
T= CMT GeEGel) XMM4 
TF OM1GeE S04) XRRZ @XRRICKIOD 

TS CNL a eE 204) KY ex4a4Icarap 

TE CMTQcEIeS) KRRICMIGISOKRAICMIQ) SAAT 
TE CMT 2 eE 203) KYMI CHI IV SeKMM TCM G) SIAL 
TF O42 GeE 206) KRRICILIISHRAICMIDIZPAT 
TF O82 a eE Seb) KMMICHTQIEKEEE CHT 54T 


CONDO KARL CAT ad 
CONLeOKaMITO4TG) 


TE COST Ade EF OP CANDO CATQe5zt 05)? KRRRSKIRICHI 2 TAO 

T= COSTA eT Gel) cANDeC41Q05E 25)? XMMMSXMMT C41 Q)2 9 TAUWO 

LF COST AU aT I el) eAND 0 C41005505)) KRARSKVWAL CHT Ge Tauad 
148 


cee wwwF (2) AND § C3) emcee ew crwemwaweceenens as eeweeces es cons cess eee nes woos 


> YIETOSKYSIAMMAS*ISCSKYCSIC/ PAMMASSIAFLAUSOCUBCLA 


LF CCSTAUCE D201) e AND 
TF (1JseEQe7) JRITE 


eo (41G0eGEe5)) XMMMEX4MIC4IG © TAUEL 


CSo71L) CHRIICII 91 3194) 


LF CL seTGe7)d WRITE C5972) CKHMAI CI) ol H199) 

Tl FIRMATCIOK a wmeeKRIICTI elLHle4 IF FC2) AND FUE) AVE ewe eg Gl E 14079 2X)? 

TZ FSRIATCLIOK eg emo K MMI CL) oleleg 4 JF FC2) AND FC3) AVEmmm ee GC E14e7y 2K)? 
KRR = KRREKRARR 
KAM = MAME KM 

331 CONTINUE 
DemwenCALCULATIIN OF HICZZTAL) erceeee, 

XSIPL = XXlele 


XSINB = XKL°¥XSCL) 
XSIMF 5 XKL-¥XSC3? 
XSIMC = Xxx1-Y¥xXS(2) 
KSI212 = XSIPL#*2 
KSIMB2 = XSIMS##2 
XSIMF2 = XSIMF eed 
MSIMC2 = XSI4C##2 

RRA = SGRTCXSIPL2¢vy12) 
RWS =S SARTCKSIMS2eVYL2) 
RRC = SGRTICXSIME 2¢yv12) 
RID = SARTEXSI4C2evr12) 


THIA = ATANCYYL/XSI71) 

LF OxSLPLecte3e) THLA = PALeTdIA 
THI3 = ATAVCYY1I/xXSI43) 

LF CXSIMBeLE ede) THIF = PAL*THI3 , 
T4I> = ATANCYYI/KXSI47) : 
TF CXSIMFeLE oe) THIC = PAL+THIC 

THI = ATANCYYI/XSI45) 

IF CXSIMCeLEeO0e) THID = PAL*TAID 

R27 = SGRTCRRACRARIPRAIC/RARO) 


THITL = eS eCTHLACTHL ae THI Ce THL OD 
COTH1 = COSCTHITID j 
SIT4i = SIVCTHITID 
F2C0 = RRL*CXRReCOTHL=XMMeSITAH1) @AL= AL H 
FECT S RRL ECKARGITALEXMM eC ITAL) oxK<x | 


IF CIJveESe7) © (2) = F260 
LF (IlJetQe7) F(3) = ©F3C9 
IF (luefGe3) GO TI 343 

LF (luefGe9) 50 TO 341 

IF (IlJeFGeld) GI TO £42 | 
IF (CIuefGell) GO TO 343 

T= CideESelZ) GO TI 344 

TF (lJeEG013) G2 79 245 | 
I? (LseFGel4) GO TD 345 i 
If (LJetael5) 59 TI 347 | 
TF (IU06&Ge463) GO TO 403 
TF (10506009804) GO TI 904 

IF (IveEGe905) GO TO 905 
I? (Lde£G0406) 52 TS 406 i 
TA2G = TANCYXS(5)65444A) ' 
KIMSS XRRICLIECITAL © KAMICLIESIT42 ! 
MSKTE KIRLOLI@SITAL © XMMECL COTA 
ROLERVLACITAL 

QS1SRA1TeSITH1 

KRSG=XRRICSD*KRRICE) 

XASGSK4MT CSP ONNIICS) 

RAHSSRCTONRSSRSLOXG5S 

R45SSRO1LMMSOORS1OXA5S : 

20205) 2 PWeTABSORIAS 

PC295) = C2 sSI/PATORRSEHOTHOA 
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203e5) = eVRLOTASSHKSKC 

PO3e95) = OCT eSI/PALOTA2ZGCRMSSHOTHDA 
B>Vs2eePAleCyerVxS(5)? 

PCAs SIERRLOCKRALCLI*COTHIMKSAL CLI ASI TNL /B2 7 e8RSGOeD TAOS 
POSeHISFRLACXRRIGL #SITHL OXIA CLIO CITHL) BP Ye24ESeOTADS 


lJ = 3 
YXSCL) = XCLIODELTC192) 
30 TO 330 
390 FP2 = F2cd 
#23 = =3C)9 
IJ = 3 
YXSCD) = XCLI“DELTC291) 
33 f2 330 
341 PO2e1) = CFP 22F2CIIS (2c *DELTC201)) 
FO3—92) = CFIZLEIWCIN Cee DEUT (2919) 
YXSC1) = x¥C1) 
YXS€2) = x1C2ZVOFELTC292) 
Iv = 10 
3G T9 336 
3492 F722 = Fac) 
F23 = F309 
YuSC2) = XCQ)*DELT(2¢2) 
TJs11 
32 TI 330 
393 2C2q2) = CHP QF ASIII C2eeIZFles2)? 
PC3e2) = CFE3=FICODS CZ e#DELTC292)) 
YxS¢2) = x2) 
YxSC3) = AEDI*PELT C293) 
Tu = 12 
329 TO 330 
334 7P2 = F2C) 
F23 2 F309 
YXSO3) = XC3II“OELTC2¢3) 
TJ = 13 
69 TO 330 


39S PC2ekd> = CFPQFACIIIG Ze MDE~T 4293)? 
P (393) CEPS=FICOIC2e# DEL TC 293)) 
YXSCH) = XCOGD*DELT C294) 
YKX503)2K03) 

TJs14 
32 TQ 330 

396 FR2=F2C0 
F23ysF3C9 
YXSC4) = XO4d-DELTI204) 

Ty = 15 
32 TQ 330 

397 PC2—94) = CEP 2=-F2CIII(C2e*DELT (294)? 
270394) = CFIZSF ICI C2 eM DELTC2049)) 
YXSC4) 2X04) 

Tus403 
YXSO7VEKCTISCOLELT C207) eaBSX7 
39 70 330 

433 F22sF 253 
FP3=F3C0 
YXSCTISNXCTISOELT(2 07) 2#ARSX7 
IU=404 
30 T3 330 

929 PLAg7ISCEPZ©FICII/ C2ee DEL T6297) 2 ARSC 7) 
PCS 7TISCFIZ“FSCOMS CZ e*# DEL TI 397) 2 ASSKT) 
CXSC7TISEKE7) 


150 


; 


133405 
YXSCBISX08) *OELT(2 08) 
$6) TO 330 
995 =>2=F2cd 
FI3sF3C9 
Tds405 
YXSCBI=SKCSI-DELT (203) 
6oO TO 330 
996 PC 29BISCFI2—-F2C6994 02 ot DEL TC 298)) 
FOSeBISCFIZS“F3CII/ C2 e*DELT(3e3)) 
yxuSC8d 2x3) 
BR emeowes (4) ce cewececeecc eres meee ewe cele: ceenerececes concesecocer = coce 
LJ=16 
YXSCLISKCL) SOZLT C4eld 
1399 CALL OFSI¥2¢ANS2) 
TF (i velGelS) 59 TD 313 
TF(IUveZQel?) 6C TD 51% 
Til JeEGel9) GO TO 575 
Le (lvelael3? 30 TI 315 
T¥(IuefGe20) 60 TI $15 
TF (idselde21)2 30 T2 317 
IF Ci veEQe22) SO TO 518 
TF (Il Jelae23) GO TI 521 
TF (I veZQe2%) GO TO 522 
TFCILvefse25) GO TI 523 
iF ¢Iv052025) SO TI 524 
Todt ded Ge261) GO TO $241 
T?€LUefGe2S2) GI TI $242 
CIuelGe407) GO TI 407 
(XseTG0808) 30 TI 868 
CT def Ge499) 59 TI 409 
(Ide EG0410)2 SO TI 410 
513 ANS? =4NS2 
1J=17 
YRSCIVSKCLI“OELT C491) 
GO TO 195 
314 ANS224NS2 
13215 
PCG LIS=CANSPHANSQII(C2Ze# DEL TC 4e1)) 
YuUSCLISNCL) 
$0 TO 195 
$75 ANSF=ans2 
F(4)S=(BISS<ANSF) 
13219 
WKSC2VEKC2) SOL TCO e 22 A3SCKC2)) 
39 TG 195 
SLE ANS? PZAaNS2 
Tu=20 
YXSCQ2ISKC2ZPSOEL TCH 92d eABSEKC2)) 
39 TO 199 
315 ANSIGQZANS2 
OMe 2dS—CANSPIMKANSIIIS CZ oP DEUTCH e2VOABSCKE2))) 
YXSC2d5x¢2) 
bJe2) 
YRSCIISKCIIOOELT COs 3d exC3) 
33 79 199 
$17 ANSIPZANS2 
TJ222 
VKSCSIERCSI“OLLT CO e Sd OKC3) 
GO TO 199 


wa ceR eee Re bee a Le cae Renee ee 


513 ANS1Q=aNS2 
PO4e SPS CANSIPHANSISI/S CZ eH OELT (4 e353) 0X03)) 
YuSC3PSxXO3) 
TJz2s 
YRSCODSXCSIOOEL TCO 4d PABSCKCG)D 
32 73 199 
$21 ANAZANS2 
tJ=24 
YXSC9I=SXC4) -DELT C494) eABSOXC8)) 
33 79 199 
S22 ANBSANS2 
POG 94) S—CANASANSI/ 020 *DELTC 499) tAZSUX0 4) D) 
YKSC4) 5X04) 


Tvs2s 
YXSC52=2005) eDELT (495) 
38 TO 199 

523 3vasavs2 
TJj=25 
YXS€S) 3X05) -DELT( 4_5) 
33°73) 139 


52% SNBSANS2 
Poe Ed EOCINASBNBIS(C2e#DELT 64953) 
YX¥SC5)=K C5) 
F542 (4) 
yx$¢5) = «¢5) 
YXSC5) 5X06) +OELT( 495) 
TJ=261 
$0 TO 199 

5291 3NAZANS2 
1J=262 
YKSCHIZSXCSI“DELT (495) 
32 19 139 

52°2 SNSSANS2 
PUG GPSHCIVAM3NSI/S (C2 eeDELT C4952) 
YXSC6)=xX(5) 
TJU=407 
YXSOTIEXCT) OOELT C497) ©ABSX7 
32 73 199 

427 ANAZANS2 
YXSOTIESXC7ISOELT( 497) *ABSX7 
T3403 
50 TO 199 

438 ANBSANS2 
PCS 7TI=SCANAPANEBI/S(C 2e#@ DEL TC 497) * ASK 7) 
YXSO7) S407) 
YKSCSIEKCZIOOE T6493) 
TJU34903 
39 TO 1399 

409 ANASANS2 
TJ5410 
YXSG33 2X03) DEL TC493) 
39 TO 199 

SLO ANB=aNS2 
PC 4—eB) S=CANA*ANED/(20#DELTOC4 eS)? 
YxS(6)=2xX09) 

eorenF (5) AND (7 eww mme ccc cewencmces wees ce mewecewews ce ecceee coeaesocoes 

THIS SUSIOUTIVE FINIS THE END PIINT JF CavITve 

IJ = 27 - 

BLS CALL CAVITY CXCENOsYCINO) 
r31l223 


r 
. 
- 
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322 


322 


824 


826 


827 


CALL SHAPECXCENOe TUPPER SBETASISII2) 
T2C1Ueb9e27)2 SO TI 323 
LE CI veEQe28) SO TO 8321 
TF (I seTGe29) GI TI $22 
LT? CLJseEGe30) GO TI 323 
LF Cl seEGe31) GO TO 324 
LF (i seEGe32) G3 TI 325 
TF (CIdeEGe33) GO TO 826 
F (IyetQ034) GO TO 827 
TP bIveEGo341) GI TO 330 
TF (L50e£0035) GO TD 328 
I? (1 JefGe35) GI TD 323 
T= ¢IJetGe37) 60 TO 340 
TFT def 20359) GO TI 341 
IF (IdeEGe4l1) GO TI 411 
TF ClseEGe412) SD TI 412 
TF (1s eEGe413) 30 TI 413 
LE (I seEGe414) 39 TD 414 
FCS) seCXCIVO“CAVLEND 
FCT) & <CYCEND*YUP PER) 
lJ = 28 
YXSCL) = XCLISOELTCS,1) 
39 TO 815 
ANP=XIENO 
ANP7TZYCENI“YUPPER 
Iv = 23 
YXS€1) = XOLIODELTCSs1) 
32 T9 915 
PC5 91) =CANP=XCEND) /€ 20 #ODELT(591)) 
ANG?TSYCENJ“YUPPER 
PC7T912 = CANPT=ANQ7)/O20#DE.TC5e1)) 
¥USCT7 Ss X02) 
YKSC23 = KC2IFOELTOCS 9232 A39S0K02)) 
Tu = 30 
39 TO 815 
ANP=xXCEND 
ANPTZYCENIOYUPPER 
YXSC2) = XO2>*DELT(5 02) #A3S0X(2)) 
IJ = 31 
39 T9 315 
POS 2)=CANP@XCEND) AC 2e#DEL TO S92) eABS(X(2))) 
ANQA7TZYCENI“YIPPER 
PU792) = CANPT=ANG TIS CQ2e#DEL TCS 92) eAS3S6K(2))) 
YxSC2) = x€2) 
Id s 32 
YXSC3) = XCZIODELT CS 3d 0x3) 
339 TD 915 
ANP=xCENO 
ANP TSYCENIMYUFPER 
YXSC32 = UOSISDELT CS sy 3) eKC 3) 
Wd = 33 
33 f2 915 
PO5_ 3) SC ANP@XCEND) /0 20 *DELT (503) 9X (3)) 
ANG7SYCEND=YUPPER 
POT%—3) = CANPTOANGTI/ C2 ee BELT (Se 3) #xC3)) 
IJ 2 34 
YXS¢3) = x3) 
YXSCO) = XCGDODEL TC 5594) #ABSCKC4)) 
33 TO 315 
ANP=XCEND 
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ANPTSYCEND HV SIPER 
YXSCAY = XCOd“PELT C504) ABS 8) ) 
TJs341 
6) TO 318 
330 COVTINJE 
FOL GPS CANPERCENID CIA DIOTC SoS eABStN6ard) 
ANGTZYCENJ“YUPPER 
720799) = CANPTOANITIAC[ePDE TS 4) eA 3SCKCA)DS 
YxS(4) = x4) 
YXSC5) = xK€5)2DELT (395) 
IJ = 35 
33 TO 315 
325 aANnesxsenvo 
ANPTZYCEND-YUPPER 
YRSUSD 5 005) eDE.T05e5) 
IJ 235 
35 $3) 918 
325 PUSeSdSCANAMKACENVD) SU 2ee CELT 6 5e8)) 
ANGTSYSENDHYUROES 
2079S) SF CAWPTOANITI“SN 2S o@ DELTC5 057) 
YuSESVSKCS)D 
YXSCGVSKCSPOOEZLT C505) 
1u237 
33°73 315 
34aC ANPSKXSEND 
ANPTSYCENIJ©YUPPER 
VESCGISKOCSIMDELT IS 95) 
IJ2368 
33 TO 815 
BOL PCSeGISCANPHXCEND) S020 DEL TCS 95) ) 
ANGTIYCINISVUPPTR 
PU7T96) = CANPTHOANTT) SC 2e#JELT( 55 6)) 
YXSU6ISK06) 
Losell 
YXSCPPZKCT) MDELT CS 97) #ABSX7 
309 TO 815 
912 ANPZ=XTEND 
ANP TSYCENI MY UIPER 
[yza12 
YXSC7ISXCTISOELT CS 97) ABSKT 
39 TO 9315 
G12 PCS e 7ISCANPXCENOI/SC2e#DEL TC 557) ABSX7) 
ANQ7ZYCEND TUPPER 
POT TIS CANPTHOANITI SC 2c DEL TCSe7) *A8SXT7? 
yuSO7)2x07) 
132613 
YXSCBISNC3I*OEL TOS o5) 
39 73 315 
913 ANP=XCEND 
ANPTSYCENDOV UPPER 
TIS9219 
YXSCBISX03) -DELTC5 93) 
33 FO 515 
BLO POSe BIZ CANP©XCEND) AE 20 *DELTC 598)? 
ANQTSYCEND@YUPPER 
PC Te BIS CAN? TOANQ7T) S02 e*DEL TE 598)? 
YxSCa)=K¢3) ‘ 


De ew ee wcel (5) cocccene 


LJ=80 
950 Caner IJFSIMSCANSS) 
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354 


3335 


333 


IF C1UeE92090) 60 TO 851 
EF CLdeE F041) SO TI 352 
TF CI veEQ042)? GO TI 353 
LF (0 Jv0b S083) GD TI 35% 
TF CI UeEQe44) GO TO 355 
T2 C1 veEG0495)? GC TO 555 
LF ClJvetae45) 50 TI 357 
Le CI VeEQe47) GO TO 358 
if C1 U0E 9048) GO TI 359 
IF GI JeiQ04%F9) 50 TO 360 
Lit JeEQe50) GO TI 351 
IF (lJefGe5l) GO TI 352 
iF (I dvetG052) 50 TO 553 
C1d0ef 20915) 52 TI 415 
CIlvefGe416) SO TI 41€ 
ClseEGe9l 7)? 30 TI 417 
CT seFSe416) SO TI 415 
F(6)s°CANS5*51G6S2) 

TJs41 
YXSCLISKOLISOZL TCS el) 

32 T2 350 

aANPZANSSE 

1y=42 
YRSCPISRCLISOOLT CS old 

39 79 350 

A oolPTCanreANS5 IAC eH CETC501)) 
YXSC1ISKC1) 

Tuses 
YXSCQ22SN(C22¢DELT (S02) 

33° T2 380 

ANPZANSS 

ljsee 
YRSCQ2VERC 22 20EL T6502) 

39 TO ase 
P(HeQPECANVIHANSSI/ C2 POEL T C5 92)) 
Luzat 

yxSc2) 2002) 
YuSCTISNCTIODEL TCS e3? 


oe cee 
naw 


33 TO 85¢ 

avrsanss 

TJ=86 

YeSUS SUC HOOLT S93) 

30 TO 3580 
F(GeSISCANPPANSS IAC 2 PIELT (593)? 
Lyse7 


yxSC30 2x03) 
TESCODVEKCADOOE TS 08) 
39 To 6590 

ANPZANSS 

luzaa 
YUSCODEKCODHOELT C3509) 
32 T2 350 

Pe be Od TC ANPHANSSIS C2 POEL T(698)) 
iPETy) 

YuSCadSx(4) 
YKSCSIANCSIOOELT (5 e5) 
32 T3 680 

AWPzAaNSS 

145230 
YESCED SHC EDVOELT(C G05) 
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33 TO 350 
851 PCG 5 =CANPHANSSI/(020*DELT(59e5)) 
Yx$¢(5)>x¢5) 
YXS€6) 2X(6) DEL T(6e5) 
; Tj252 
33 £9 850 
852 ANPZANSS5 
¥XSO622X(5) -DELT(S 95) 
lusS2 
32 TO 850 
353 PC be GS = CANPHANSS) / (C2 e*DELT(S95)) 
YXS¢6)25x¢5) 
TRSCTISXC7ISOELT C597) #ABSK7 
12415 
33 TO 350 
415 ANPZANSS 
YXSCTVIESXCTISOILT (607) @ASSK7 
TJ2615 
32 TO 359 
B16 PCG TIS CANIMANSS 4 C2 o# ELT C507) A387) 
YXSC72 5X7) 
YXSCBISXCSI*DELT (S93) 
tus417 
33 T) 880 
$17 ANPSANSS 
YXSCBPZXC3)$OELT( 593) 
135413 
Go To ase 
913 PCEeBISCANFAANSSI/ (2 e*DELT(598)) 
YXS(BI=X(83 
Se ewnl (9) ewmmameece 
SIDZSINCDELTAD 
COOSCISCDELTA) 
SSASVKS (GO) #STIOVXS (7) 
CSB2YxS(4)*CO0 
, SSASYKS 64) *STD-¥xS (3) 
OSB82YxS$ (4) "C09 
I SSAZS°SAee2 
CSB2=c56e¢2 
ISASsISASe? 
ASB25IShee2 
PILSYNS (4) COD—-2 ee YXS(73*eSID*#C09 
997=2C5A2CS832 
PUQSYXS (C8) ero d—-2 oe YXSCB)@SII9eC99 
9322395429832 
VAISS aITCCCSaV]|eC$32)/09S8SA2¢9532)) 
ACO=ALOGC(RCERD) 
SITO=aATANCZSB842SA) 
SITS SaATAN(9SB5/0SA) 
TE CSITOeiL Sede) SITIZ7ALSSIFO 
Scossitcesito 
FCB S—CCOPHMACISSIIeSISDI/PAL MSINCALFAL*©YXSC5) DS COSCYXS65) GAMMA) 
70891206 
PCBs 2) =00 
P(Bs3)206 
PCB SP=0PJ1/ P81 =PU277292)/0aT 
SUBe SPSCISCALFALOGAS4AD/STIS UY KS C8) OGAMHA) 202 
PC3 062206 
FCB TISYNSCTIMCOSCOILTAI/SC PALO COSA2*CS82)) 
CBee BISeVKSCAVMCISCIZ.TAISC FAT (0542 09852)) 


ee a ee 
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NCAVISNCAVelL 
33 253 ICV=leNCAVIe2 

253 wRITEC60252) CAVXCIZTVIeCAVYCICY) 

2352 FIRYATCLOX es @CAVK Se eF LO eS eo SX aM CAV Y=eeF 1065) 
33 123 ITx=1¢3 

129 dRITECSol3512 LTXeF CITX) 

L351L FIRGATOAUXe24F Colle lid tel lee?) 
09 132 IUP=108 

132 MSVITERS oL3S32 TUPg CPCI S2eJU2) ep UIPH198) 

L33 FORMATCIXe2HPCeol le Ste JI=eSCEL3 e591X)> 

335 CONTINJGE 
CALL JETERMC Pe Se DE TIO) 
20 25 ITOET=I199 ! 
33 26 ~?G2lea 
GQOLOEGe IDETISEPCLPGe IDET) 

2$ PCLOSeLIETI=ZF (LPG) 
CALL DETER P e be DETID 
TFCIDET eEQel) DELZ=HIITE/IET39 
TPCTOLTeEGe2) OLLIZIETESOSI 39 
TF CIDET eEue3) OELDAIZTE/DETIO 
TFCLOET eEi04) JELESIETE/IZT3I 
TPCIOET CES eS) OELFAITTE/OETSO 
FCIDET eESeS) DELSZIZETE/DET3C 
IS CIDET ef de?) OLLASIETE/IETIO 
TICIBET eEGe9) DELI ZIETE/OST39 
33 27 LPGs198 

27 PCLPGelDET)=GCLPGelLIET) 

25 CONTINJE 
AGL) SXC12¢D0ELB 
KCZIEXC2I4D ELS 
KOS SKC32435L5 
RCO SECAIFDELE 
KC5)V 5005) +550 
XCG6I=K(HBIODELG 
XO7I SKC 7IOIELA 
xOS8YSXCAI*DELI 
39 50 SYN=198 

$0 WRITECS e961) LYNG XCLYW) 

SL FORMATCL Xe 24XCotle2s) 2961407) 
AGSSzARSCIZOB/ X10) 
AZSCHASSCIELC/XC2)) 
AZBSISASRSCIELO/X(3)) 
ABSEZABSCOELE/X(4)) 
ARSF2ARS(IELF/XC5)) 
ASSSZARSCITIG/XC5)) 
ASSHSARSCIELH/XC7)) 
AZSTSABSCIELI“SK(3)) 
<210=0 
TS CAbSBeLTeSTILY KELIE1 
TFCABSCeGTeSTOL) KEII=0 
TFCA3S5 eGTeSTOL) <EII=0 
TFCABSE eGTeSTIL) <£I920 
IF CABSS eGTeSTOL) KEIISO 
Li CAGSGeGTSTIL <EII=9 
TF CABSMeGTeSTOL) KELO=0 
I CaBST eS eSTIL) <ELI=20 
IFCKETOeE2e1) 69 TI 35 
TeIe1 
WVUITECSe42) JT 

92 FIRMATCZOX9 1 GHITER ATION NOotel 2) 


a eS cca ie 


ae ee ee 


tr7a) 


TSCLelae4) 39 TD 35 

$0 Tu 55 

TF (leZlaeM) 50 TD 36 

GO TO 33 

WRITE( 5937) 
FIRMA CL Xe SAHIXENE a SID NOT CONVERSE a@ITHIN t41T) 
TF 0X07) eSTeSI7) XC7)2SI7 
BRECKCLe&LTeSIL) KCLI=ASTL 
ST102X€1)¢20806 

TF Ox 62) ok FeSIT10) (2225710 
SI112*K€ 22920006 
TEOnOS) ee PeSI11) xXC3SST1I1 
SISAXCZ eZee Ol wT C3 eld 
LTFOXCSd cL TeSIB) XCBI=2ST8 
LECX(Sd elielele3) X05)21 05 <3 
TFOxC4) eLTeSIS) XxX¢4)=SI6 
SIS=2CeS OAL =-GAMMA) el Le 002) 
1?€x€5) el Tee) 63 TI 51 
T5Ox5)2 e57TeSI15) x05) 2815 


33 TQ 382 
TF CABSOX(5)) eGTeSI5) x€5)22-S15 
CONTINUE 
YEQ2LLeSCETIRG PeOeELod7(XC FI 
) wwmmmmoemmene REZ NAAT SSZL EMACEB P70 Xr emer ew nme me X20 TAMROF 2221 
RETURYV 


ZNO 


ore os | 


30 


ae eter pane or enet 


SISVOITING JFSIMLCANS oNOF eXTA) 

OIMENSION xST(8) 

COMMON YCOCeSBETA2 

SIMMOW XITHC2002 eo XITVC 200) oe ANSG2SC239) oe SARC2¢209) 

CIMMIN CAVXCLOOSCAVYC LOO se SETASSZETAC eo KCCCeNZAV oLPMMeNS2 

COMMON AGOLIO) oTSHAV2 gNCAZV 9 FSTANC1IO) pSBTAN26199) gpSETANZ 06106) ‘ 
SOMMON FLASANG DELTAS ISAP ALF AL CSAMMA 

SOMYON SSETA oKK4e I COL ey SARTOIC513) 

COMMON IDUL eo XAaXBeXCeTANGoE MoV Co Ve JBIGS eXL315Se B81 GSeSMALS OSS 
COMMON XSWCB De See ERC OV VV eXMe ITER AG SKXSLOC A) eSxXSII0CS8) o YXS (8) 
SIMMON PSI ZeLPeSARCC 513) eSAICICS ISI LP 4e OE 

SIMMON BETANC51399 BI TAMC513) pI Ue lLPKe XTIC 200) 9X UIC 200) 9X0x 

CIMHON KRIIVOGAZAA 9323390252 

COMMON AAAAe BBBBeCOIle AB edhe C Bed Be TSAUSC(100) of GAUS (100) 9NGAUS 
CIMMON/SPCJESPACE 


NOF = C CALLED FRIM FLINTe 

NIE = 1 CALLED FRI WINNT FIR VEAL FART. 

NOF = 2 CALLED FIN4 TINT FIR L4AGe PART eo 
NOE = 3 CALLED =294 CAVITY OXSNIW AT FC5) 


TF Cll>Lertaed) 59 TI 3 
30 10 [@ = 198 

XSTCIS) = XxSNCTQ) 

3) TO 12 

33 11 TH = 108 

XSTCi4d) = xSel4) 


CONTINUE 

LECITIRAet2e1) SO TI 222 

60 TO 223 

32 22% IL< = low? 

FOTANCIL<) = SBETA 

CONTIVJE 

CSPACE = CLo+KSTCLIISFTLOAT CIP KM) ~ 
FSPACE = CSPACE/FLOATCUPM=L>OK) 

oF MBS P4O3 


XBET = -le*CSPACE*FLIATI(LPK=1) 

XSTiselerCSPace 

321 3 35TAvC2) 

B21 = OXSLLSKSTCQAII/CCXSTLeLed eC XSTCOLIMOXSI1) #OXST12xSTC3))) 


8715 = SQzrcaP1) 
F3 = 3£1ea71S 
X01 XSTCOIOSINGIELTA? 


nants casts ae nate a nineties ante natin Bnanisaelineah annie eseinine aaiett tte ne 


YY1 = KSTC4) *COSC(DELTA? 
vrs 5 YY¥lee2 

PLY = XSIL -xKX1 

PwM2 PL¥ee2 


PoMA = PLA2QEYTI12 a 
>xSR = PLY/PLMA 

PxSI = YYL/>LMA 

IFCNOFeEGe1) F3 = F3e>XSR 

LE(NOFeEGe2) F3 = F3*2xXST 


TFONOF ofS03)2 SSSFICCKSTLMxKCA) 


ANSA=06 

39 1 L = 29 .PM%3e2 
F, = 73 

S?ASE = CSPACE 


IF (1eSfeLPX) GI TD 30 


XSI2 2 ole tSPACTeFLIATCL) 
KSE3 = xST2¢S9aCzE 

$6 TO 31 

SPATE = FSFALE 
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ve 


$i 


KSI2 = KBLT+EPACE*FLIATCIM“LPK+)) 
XSI3 = xSLl2eSPact 


322 = SETANCI+1) 

3Z3 = ZETAVIIO2) 

A272 S CKSLASKSTC2) DA CCKST 26h od MC XSTOL) MKSI2) CXS L22KST 0399) 
APS = CXSI3S“XSTC2)2/CCKST 341 0) @OXSTCL) @XST3) eOXST3-xST(399) 
A228 = SARTCAP2> 

AP3S = SQRTCAP3) 

F2 = 23E2°aP28 

FR = 3E3*aP3S 


WA2 = XST2-KXL 
MA22 = HAQee2 
HB = HA22¢YY12 
HOR2 = HA2/AB 
4C12 = YYI/HB 
AAS = XSI3“XxX1 
mA32 > HASee? 
HIZSHAS2¢YY12 
NSZZ = 4A3/ 40 
acl3 = yvi/4o 
TE ONOF eEQel) F2 = F2e4CR2 


TFCNOFe@E Qed) F3 = = 3eHCRS 

LF QNOF CES023 FZ = F2eHTT2 
TFONDFeEie2? FR = F3e4cT3 

TF CNOFSEQ6 33 F2 = F2/EKST2=KCA? 
LEQNDF 6803) FS = FSV OKST3SxXTAD 


FSUN = (CFlederF2oF 3) eSPACT/ 36 

ANSA = ANSAFSUS 

CINTINUE 

SIL = SURTC CHL ee XSTIC 202/61 ee XSTCS93) 

S22 = SUATUKST EL ei) 

SG3 = SQRICCXSTELIHAKRSTC2I/OKXSTCLIOKSTO3))) 
ANTI 2 ZETANCL I e2e eS 2ITECSPACIeSai/sa2 
ANT2 = GETANCLPMD® 2eeSQRTCFSPACE #S23/SE2 
APLA = *1eeKX1 

APLAZ = APLAw#2 

473 = XST(1I=xXXI 

APL52 = Ar LZe€¢2 

IF CNOF eEGol) ANTI = ANTLEAPLASCAPLA24+YY12) 
BF ONDE eEGe2 2 ANTL = ANTLOEYYL/CAPLA2+YY12) 
ITF CNOF eEDe1) ANTZ ANT2*APL3/CAPLI2¢71123 
I? (NQF eESe2) ANT2 ANT2eYVL/SCAPL32¢°YV12) 
TF GNQF& eEGe3) ANTL ANTI# Ce LesKCAD 

LE ONDF oF O03) ANTS = ANT20CKSTCLIE XT AD 

ANS = AaNSAeaNTLieant2 

RETURN 

EVO 


e oe fee 
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SJBRIIJTINE JFSIM2CANS2) 
DIMENSION XC3 eXITOSP OV VCSIOXITICS I oEXUC3I oF CNSC3) KX ST(3) 
COMMON YCCCeSBETAS 
SIMMON KITMC250) eX ITNC200) oANSS2S06290) oSARC20200) 
COMMON CAVXCLOO eCAVYC 100) sdE TAS eSETAC eXCCCeNCAVeLPMMeNS2 
CIMMIN AUOLOO) oI SAA ANEATVe SSTANCLIO) SBBTANZC190) oeBETANZ(109) 
COMMON FLAPANeDELTAs OSAP eALFALSGAS4A 
ca4qoy SSITAgK «4g ICPY 95ART699(513) 
COMMON TOU Le XkoXSe Xl eo TANG EP oY CoV %oeIBIGS oXLSITSSe BIESoSMALS oDSS 
COMAON XSN(5)9 SeDoERC GV VV oXMeITERAeSXSIOC 8) eSXSIOOCE) o YXS (8) 
SIMMON PSLZe LA eSAICES13) pSAVCICSI3>»  PMe DE 
COMMON BETANCSL3SVoBETAMCS13) oLve LP Ke XIITC 200) oXdI(20090x0X 
SIMAMOV XRIINDGAQA4 93233 eC20° 
COMMON AAAAs BBBBeCClOCe AB eS Be lB e DS eTSAUSC(1900) peGAUS (100 IeNGAUS 
SIMMONS SPOSISPACE 
32.13 1651404 

L3 XSTCIS2 =¥xS¢16) 
741=301415392653 
CCCISALOG CL e+KXSTC(6I/C 2eePAl) 
JSIZECISCALFALSGAMMAI/SOSCKST C52 ¢GA44A) SESPATE 
KAKKSALOGCUI2) 
SCSPACESCle® XSTCLIVAFLOATOCLPN) 
ACSPATSOeSHCSAACE 
TSPACESCSPACESFLOATC IP MML PK) 
AFSIALTS0 oS PF SPACE 
KSET=H-Le+CSPACE*FLOATCLPK=1)3 
CIESCISCIELTAd 
SJESSINCOZITA)D 
SASXSTCL)©XSTC4) #SDE 
S32SXSTC4IHSIE 
P2PSCOE/(GAS S245 Be #2) 
FINS CSS SDSAP MPP IMXSTCLISCPAT eSOWTC1eeXST(6))2 
LPKISLPM=LPKey 
33 1 LP=1_9.7M 
TECIPeEGel) GI TO 2 
HSPACISAFSPAC 
SPATESS SPase 
TECIPeSTe.F KI? SO TD 30 
XCLISRSTCLIMSPACECFLIATCIA2©2) 
KC2EXCLASPACE 
KO3VSKCLISSPATE 
39 TO 31 

30 HSPACE=HCSPAC 
SPALESCSP ASE 
XO1) SXBET@SPACE#FLOATC IPeLoKI-1) 
XC2) SKCLI9A4SPACE 
XC3ISKCLI SPACE 

31 FINS CLISFCNSC3) 
N33 
TFCIPeEG@eL PM) NK=2 
33 9 Le2ev< 
TF yeEG0407) GO TO 7 
TF CI vel 20408) GO TO 7 
ITF CT JeEG0409) GO TO 7 
Tit selGet10) GI TI 7 
TFC JeGEe23) 59 TD 5 
GO TO 7 

3 TFL el de2) KITC2EXITAYCLOMM[ 261) 
TRC oE Ge 3) KITCSIEXITVCLPMeT P01) 
$2 TO § 
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el tt atm atessaeeninainsaaacsssacsaamaiaaiaiciesinan ttl 


a 
- 
- 
- 


T 


351 


CONTINUE 

YvCT 2x Cr) 

43 CA.CJLATE Gl 

CALL OFSIMSCYYCIDe XITSCI ele) 
MITCIPSxXITCCI) 

LF CI UeEGe18) GO TO 5 

63 TO § 

T2CLeEGe2) XITMCLAIM@“ITPOLISKITCT) 
T7(relGe3) XITNGILPMeLFeL ZKXITCI) 
CONTINJE 

EXUCTDSEXPCHXITCI2)D 
Sf2xCI-xSTC4) 2SOF 

GO=xXST(49)*COE 

PXASGl ee 2¢5) 222 
DHONSIGAPHXCTIMCOE/C PKAePAT? 
FON CIDP SEKUCI) *OwIX/ S52 
L7Cn CL ebloded FCUSCLIS=FCNSCLD 


8 CONTINUE 
E2< If FONSCID IS ALMAYS POSITIVEe 


2 


1 


IF CLP efGe0.P4) GI TS 20 

39 TO 21 

PPRECIES COL eH XSTOG) SOE) we 2eCKST(G) COE e#2) 
SS 3S0SAP tPPQSPAL 

FONSC3) FES 

SUMS CFON3C1) *FON3C2) #4 e+FON303)) 84S7ACES 30 
ANS2SANS2¢S JM 

LEC ds e£G0l8) SARCOLIVHIP SL) =ZANS2 

32 TO 1 

SARS Cu PMI S060 

ANS2=06 

CONTINJE 


XITNCLPM)=G1 AT POINT Be 


XONT 


€1)2G1 AT PDINT x=le 

XIE TNCLIMDSSTO Le RKA KS PAT 
XITVCL) 200 

ZETIA 

ENO 


qpararvereas 


13 
12 


SJBROUTINE OFSIM3SCYe XXII eI Pel) 

DIMENSLIN XSTCB) = A0290) 

COMMON YCCCe SBETA2 

CIMMON XITMC 200) eX IT VC 200) pANSE2S( 290) gSARC2 (200) 

CIMMON CAVXC1OO) eCAV YC LINO SIETABSIETACe XCCCeNIAVeLPMMeNS2 
SIMMON AFCLOO) oISA4A2 2 oNCHEV 9 B39 TANC1I0) oBETANZ¢100) s3ETANZ(100) 
SOMMON FLAFANG DELTAS JSAP SAL=Al eg SANNA 

CORYON SBETAeXX4eICOls SARCOO(5S13) 

CIMFMIN LOD oe XA eNBe XS oTANG ep AoC a VA GISIGS oXL3ISSe 3ST GS eSMALSoDSS 
COMMON XSNCB)o CLESERCOVYV oXMeITERAGSXSIOC8) ¢SKXSIOT CE) o YXS (8) 
SIMMON PSIZeLPeSARCC513) oSARCICS13 Do LOMeDE 

SIMMOW BETANC5S13) e525 TANC513) ol VoL PKe XIIC 200) eX UIC200) oKIK 
SIMMIN MRIJND9A2ZAK03233eC20C8 

CIMMON AAAAG BIBS eCCI Se AB a S5e CZ oIBeTSASS(100) ed SASS CLOO eS NGAUS 
CIMMON/SPC/ESPASE 

CIMMIN/TAJI/SUTAS 

COMMON/CVTYL/CAVLENe SIGS2 

FUR INTESRALS TO SEF ZVALUATED 3EFDRE xI IS O3STAINEDe 
NITE THAT FRIVIDUSLY INLY OVE SINGILAR INTEGIUAL @ AS 
SALSULATED IN GCASCAD AND CASTAIEs 

Sic THE NITE OF TS 3351 FIV SIIR INFEGRALS» OIT JF ARICA 
Jad ARE OF SINGULAR TYPEe 

IT=CIC21e£&200) GI TI 9 

00 12 ITSl=148 

XSTCISI)=SXSNCIST) 

33 T2 12 

02 13 JTJU=1e8 

KSTCuJI=VKSCUTI) 

PAT=30141592653 

COCLAAL IG CL eo @XSTCE)/C2e MOAT) 

SSG=SGRTCLeXSTCE)) 
U2SCISCALFALSSAMMAISSISCKXSTCS) #SA4MAD ESPACE 

J22zU2ee2 

J22THlLe/U22 

TSO=SSG*CAVLENS IJ 2 CL esCAVLEN) 

TAUMOSALOS( TAO) 

Ta LSSAWT (Le mCLeHKST(6) -J221 eC AVLEN) 

TAUWI=ALOG(TY1) 

LF CSTAUeT20U) TAUw=TAUNO 

LF GSP AsSeD adel) TAJasTAUd) 

FIRS Jlecesese 

TF CITERAcEIe1) GOI FI $49 

30 TO 61 

CINTINJE 

39 52 TZU = 1lsLP™ 

ZETANCIZU) = SBETA 

JETAMCI ZU) = S3=5TA 

CONTINUE 

CINTIVJE 

CSPACE=Cle* XSTC1)/FLOATCLOK) 

HCSIATSOeS#SSPACE 

FSPACESCSPATE/SFLOATC.FMOL PK) 

M= SPAT SC eS5*e= SPACE 

MALT SL etc SPaCl er Jari P<el) 

AB2=SAIRTCKXSTC1> 416) 

ADZSESAIRT CCL oo VP eCxSHlLI<v)) 

A3G =.SGRTICCXSTC3SMVISCXSTI20°7)) 

A33 = A330A36 

TI2SLPM nT Pod 

busel 
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um aaa 


Li CLelde3) [USSLPM Hf P41 

T3CTeEa 00) LUSSIP 

SECSBETANCIJ3) 

TFL eEGe2) JZECHBETARCIU2) 

FAASS=IC/A33 

L241s.2 4-1 

00 1 Iw=2ebPM1 

SP ASESCSPace 

T&(CleeGTelLPK) 69 TO $5 

XSKseLe *S2aACE*FLOAT (I d=1) 

$2 TO 46 

SPACESFSPATSE 

KSKSXBET® SPACE ePLOaArC Iw 2<) 

TF CL eEGe2) GO TO 6 

Ti(ImeldelJ3) GI TI 1 
EFSSSQRTC Cle XSK) eC XSTC1L) =xXSK)) 

FSAL = SGITCaXST CS -CSKI/CKSTO2) -KS4 9) 
FS = FS*FSAl 

FACT AISCZLETANCT PS ESMeR AAC XSKOY) 
CONTINJE 

I= ¢LeEQe2)? GO TO 30 

KPLZ“LetHCSPAS 

XP 2Z=XOL*CSOACE 

X2G4=XST C1) -AFSPAC 

XP SSXPeeFSPACE 

SSLSBETAMCLIS SORTCC Le XPL) We CXSTCL> XP 1d) 
FSQ2SSETAMC 204 SART CCL eo eX72) eC KST CL) 9K72)) 
FSSSBETAMCLPM=2)/ SARTO C1Le#XPS2 HCXST CL) —XP 3)? 
CSGtBITAMCLPMMLI/ S2ATCC Lex 24 de CXST C1 I=xXP4)) 
FSAL = SGRICOXSTC2<KP1L/CXSTC3°xXP1)? 
EFSA2 = SQRTCOKST C2 KF 22/ CK ST C3) 2927) 
TSASSSARTCCKST C20 2K2 39/7 0CKXSTO3) -K73)) 
FSAGSSORTOCXSTC 22 © X2 G2 /0XSTCS5) @XP4)) 
FSISFS1leFSAl 

FS2sFS2eFSA2 

FS3SFS3eFSA3 

PSQ=FS4e=SAg 

FALSCTSLOFAAI/ CKA1LOT) 
FIQSCFS22FAAS KP 27) 
FOSSCESIMFAAD/ (XP Zev) 

FHS SGmm AAI CK 247) 

TECLUSeEQGe2) GO TO 21 

TFCL IS eF20eLFM1) SO TI 22 
TFCTUSeE Ge. PK) 52 TI 351 

FACTU3S) SOe5e (FACT US 91) oF ACI S3+1)) 

33. T9 30 

BETOS2e *BETANCLOK) MSZ TANC LPM 01) 
KJASXZETOFSP ATE 
FPWOBETOSSART CCL e@XIA) OC XSTO 1) =XOA)) 
FOwA = SQRICCXSTC2ZV9 KDA OKXST CS) 9K 08)) 
TPesFOUA eT PY 

FLP<CSCEPaeFAAD/CXOA-Y) 

FACIUS) S005 * CFACITI3+1) OFL IK) 

39 TO 30 

FACI SSP SCF IL OF P2426 

GO Tu 30 

FACLUSISCZPSOEP4D/ 2, 

X1TS06 

LPM3=LPMe3 

SPAT ESC SPACE 


15 


= Fe ee ee i te fei 


Q2 15 JA=2eL P4302 

Le tshoctei 24) SPACES SPATz 

MISMO CF ACSA Ob eeF AC JAS1LIOF AC JAZ) DOSPACE/3e 
*(Tlelde2) 30 TD 35 

ALSSS0 oS HHO SPACHC FPL OFAC 23) OCEACL AMS LD OFPR) HD eS e4FSPAC 

XKTs4le 

4J=39 

LOMASLPM<§ 

LECTUS eGEeL PMA) XCITZ2010 

TECIUS eGEeL PMA) KUZ1II 

JIZSCITTAMCLI<“BETANCLIII/XAI 
FIVSCIOTANCLPMIMKBETAMLIAMDIIALKT 

ME FOHSSPAS/XKI 

AF ASATSPACT/XKKI 

ET3SFP1 

SJ3SFI4 

X14206 

X1L1=206 

32 202 [F4sletUe2 

ETISFT3 

FJLZF JIS 
XM2EXSTCLICHFSPACOHFFOeFLOATCITM) 
MMVISKAQOHES 

XT22-Le MHSSPACHF4e= JATCIT) 
XT3SXT2en=e4 
JETAQSZETAVMCL IMT *BIVeKLOATCITS) 
SETASSBETA2*¢60Y 
JETT2Z=BETAC1) -BOZe= JATCITY) 

SE TT3SSBETT2-802 
FSQSGZTAS/SART CO Le OK 42) HO CKST C1) =xX42)) 
FSSSBETAS/SART CC 1Le eK43) MO CKXSTCLI“X43)) 
FVASBLTT2¢SSRT CC LeoxT22e(XSECLI-XT2)) 
FV BSSSTTI SSA CO Lee XPS CKSTCLINXKT3)) 


FS2A = SQRTCOXST (222442) /0XSTC3) =XM2)) 
FSIA = SQRTCOXSTC22-KASIAOKST OS -K435)) 
FV2A = SGRICCXSTC2)-XT29/¢XSTC3) =x T2)) 
FY3A = SGRTCOXST C22 KT3d/OKSTC3) ox T5)) 
FS2 = ©S2eFS2A 
FSS = FS3eFS3A 
F¥2 = Fv¥2aeFV2a 
Fv3 = FV3eFVSA 


FIQSCFSQSF AAS CKMQ7) 

FISSCFSSSF AAA UXMS HY) 

FTQ2SCFV2eF Aad CXT2-7) 

FT SSC V SOF AAS UKT 3 <1) 
XIGEXTGOHE Fe CF ULCEU2 tO etFUS) 36 

KLLSAXL Lod 4a CFT LOE TI em esl 3s 36 
XAS=BETANCLPM) #2 otSanT CHEF) /CAB2eCXSTC1)0-Y)) 
MAGA = SORTOCKST C2 RST CLIO OKST C329 KSTO1) DD 
X44 = XAG*XAGA 

XTS=xr4exas 

KALSBETANCLI e2eeSIVes ad /CRZ2 ECHL oH V)) 

KALA 2 SGRICOXSTC2) Led /CXSTC3I%10)) 

MAL = XAL@MALA 

XITTSXI1¢XAal 

MTSOXT+#KT23¢xXT 1 0x1 4) eABS/PAT 
KESXTOBECAALOG CO XSTILISVOAFOI SC Let YoHFH) D/PAT 
KAXT1S°xT 

32 TO 35 


35 XR1=Z-1let0eS*HCSPAC 
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RUVSSKVLEMCSFIAS 

AR@EXSTCLI Ge S#HFSPAS 

KVUZSKWHEMAFSPAC 

CTLEVe SCALE TANCLISCSETAMCLIIG SQITCOL eo eKU1) OCHS TOL) e219) 
ST2OeS MC ZETAMCLIOBITANCQ2I9/7 SERTOCL eo XR2) eC XSTC1) oKR2)) 
FRESH SO CSETANC AMMAR POSE TAMCL AML SQRTC CL eOK23) COKSTOLIOENRS DD 
FTG=D eS eC BETAMCLPM=1I*BETANCLIAM) ZS SARTCCL eH XXV4) OCXSTO1L) Oo XRS)) 
TTA = SAITOCKSTC[P-CALISCKSTOS) -KRTDD 

FT2A = SGRICCXSTC2d2KR22/0XSTC3) -XR2)) 

ETSa = SQATCOXSTC2d-aRSVSCKSTC 3D -KRSDD 

ETQA = SQRTCCMSTC22 CUA /CKST OS RVG) 

eT ETLe&TIA 

FT2 FTQeFT2a 

ET3 FISeFT3A 

FT4 CTaerTaa 

FRLIDCKT LMF AAS OXRIMY) 

FRISCH TZ<F AAS KR 227) 

FASSCOTIMFAAI/CKRS 1) 

ERGZCT TGF AAIICXRG =) 

MEAL SP eS OAT SPAC HCE IL OF RZD OO eS ONFSIAT CLER GOR V4) 
XIP220 e25eHCSPAC HC FR26FAC2)) M0 e25eH= SPAC HCFA CL PY) oFR3) 
XT23=Sx 191 ¢x1P2 

X“T=216 

X4IT2=426 

4Jj221 

M2s4U=2 

We FMASL FM M5 

IFCIUZeGEeL PMA) XMI=101le 

lL? ClU2eGEeuP MA) KHI222026 

TF(TU2eGEeL PMA) 4UZ1LO1 

TF Cl u2eGEe PMA M2S4J92 
ZETVSCZETAVCL PMD MSCT ACL 42D KMT2 

BESSHO eS ACZETANCLPM) OBETAMCLPM<1 >) 
4S?5=0 eS edt SPACIKAT 

FQ3=FR4 

ZETYLSCSETANCL MBE TANCLIIANNI2 

JESSISO oSHCZETAMCLI+IOTANCI)? 

HSPS1=D eS#4CSPAC/XMI 

FASLSFaL 

XT 12Q6¢ 

X142G0 

20° 129 [TLeleM2e2 

FRZ1LSF23 

FQ112FQ31 

XQSKSTCLICAHSPECFLIATICMUHIL) 

K3SK2¢hS75 

KR21s°LeeHSPS1I#FLOAT( SUMIL) 

KSLEX2 LenS 761 

BETA2=BESS*+B3ETY*FLOATCIL) 

SL TASSRESSOSETY*FLOATCIL 1) 
SITAZLISSESSL“SETVIe= ATOLL) 
BETASISSETAZI-“BETY1 

FJ2ZLSSETA2L/ SARTO Loe oXK21 eC KSTC1I9K21)) 
FUSISZETASIZ SORTCCLe + X312 eC XSTC1I 9X31) 

FJZIA = SIRT COXSTC22 9K 219 70KST C32 9x21)) 

FISIA = SARTCCKSTC2) -K3S1T970KSTCS29xX31)) 

FJ21 = FuU2L9FV2IA 

FJ31 2 FUSLOF ISLA 

FAZLSCFUZL HF AAD/ UX 217) 
SCSSLSCF ISL OF AAS EX3L OY) 


FU2ZSBETA2# SQGRTCC10%X2) #CXSTC1)-x2)) 
FJSSBETAS/ SGRTCOL eee SDOCKST CLI“ K3)2 
FY2A = SQRTCCXSTC2)2-K22/0KSTC32 -x2)) 
FISKA = SQRTCOXSTC2<K30/OXSTOC3) 2X39) 
FJ2 = FJ2eFJ2a 
FUS = FUS*FUSA 
FQQSCFUZ=FAAD/SCK2-Y) 
FOSSCFUS$FAAD/CX3-Y) 
MTLSEKCLoASPSLeCEIL Le a2iedorF asi rsde 
129 KIGSKL4edSP5eCFOLOGe eK Q2eFQ5)/ 36 
XTAZ2e*SGRICHSPS) @ BE TANCLOIMD CCAS Se CKSTC Le)? 
MIAR S SAITOEKSTCAd“KSTCLIW/CKST E32 -KSTC1) dD 
XIA = XIA®XIAA 
MIGSXITGoxla 
XITSS2ee SGRICHSPG61) eFETANCTI/SCAB2e (HL eH) ) 
XT3A = SOITCEXSTC2d.¢1L97 OxSTC32 +12e)) 
XI8 = KIBeXKIBA 
ATLEKTL «x18 
KRISOKLOXITL OKI 23¢N16) ARS />zaT 
ALEXI*BECMALOGCOXSTOC1L) eVedS2S) fC Le t¥oHSP61I PAT 
XKL1LS¢xKI1 
36 CONTINUE 
eecwene![ 2 cow cece core 
ewrmceomownT© y [TS .E58 TIAN ZE%Ie THIS IS A 
eromeAEGULAR INTEGRALe WHILE Y o5fe OF THIS IS A 
wwoerSINGULAR INTESR ALS 
BJT THIS IS TREATED AS A SINGJLAR INTEGRAL ANYSAY 


greverarnas 


Isicss 

XOASY 

Shwe TISCSRaeSMoXTASISIC) 
412252 

ARGLECXSTOCL) “YY 


TE CARS L eaT oOo) ARS THARGL 
KXTZEXXT2eAZ SAL CSUARSL) 
XXT22eKx1l2 

ewvrwe! Jeweoeosoce 
JSE CHEZSYS4HEVeGAUSS GZIAQRATIREe 
ICE) ARE ALREADY CA SILATED IN SJBVOUTINE FLINTL 
ANO PASSED ONTO HERE ZY COMMON STATZIMENTS 
XXI3 = Oe 
BPCS = (XSTC1LI*XSTC2))0 e065 
CMB5 = CXSTC2KST CLD 005 

A351 = CSPT5¢1e2/04S5 

A32 2 CHBISSeXSTISII THES 

33°1290 [834 = lsVCd3r 

MAL = LemshUCISUM) 

4h2 = CAUCISUTD CASI WCAS2~-AICISIM4)) 


eperares 


SHA2 = SQRTCHA2) 
7313 2 dAL/SHAZ 
CSAI3S = CHS5S*AUCISUMD OB>C5er 
220 XXI3 f XXTS¢FSIS“ZE3AT3 
MAKT3 S MMLSOPAT/SNE AY 
UU22 = COSCALFAL*GAMMAI/COSEXSTCS) 3 AMMADSESPACE 
: VICTSK/PAMMAGOISCTSXCSIC/IASMAGSIAF_ACSOC = 22U5 


WX3 2 CCCL©ALOGCUU2Z2)/PAT 
MKI3 2 MNT S*ABSORKS 
easeel §oowseacace 
JSE CHESYSHEVeGAUSS QUVADRATIRE FORMILASnwe 
meme eR BT TAN2¢C1) AVE ALVIZAIY SALSIJLATED IN 
SUBROUTINE FLINTL AND PASSED ONTO HERE BY 


Crpererea 
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> CIMMIN STATEMENTS 
FPCS = (XSTC3)exXSTC2)) #05 
EMCS = OxSTC32<KSTC2) 0 005 
AGL = CEPTS eLeI/FNCS 
A42 = CFPCS<KSTCLII/ENCS 
xK14 = Oo 
90 130 ISUM = leNCH3Y 
RAX = COBTANACISUMD PAT) @CLesAUCISU4)?) 
R3x = CAUCISUGD AGL) *CAUCI SUM) 0442) 
SR3xX = SGRTCREX) 


ROK SF RAX/SVBK 
ROX = FMCS*AUCISUMN)D += 2CS5<¥Y 
133 XxKI4 XXTQ @ ROXIE 


XN1S = XXI4@PAL/NCHBY 
XXI@ = @xXxT¢eaB3/7al 
DemeeelSecrccoccece 
SPL4SCXSTC7I 01 0) #05 
SYLADSCKSTC TIL ed #05 
ASLZCSM LM XSTC2) 04 C25 P 1H) 
AS22054 noxSTC1) 2/0521) 
ASS=CGMIH@XSTC3)09/0-4571H) 
XxXI5506 
33 76 TSIJMS1e NCHS? 
MAL seCle*AUCISU4)) #CAUCTSI4) +AS1) 
482 SCAUCISU4) ©A52) #CAUCISI40 6453) 
i FSASSORTCHAIL/AAQ) 
t F5a5 CoS lade aU CIS IM eS4iddey 
70 XMISSXXIS*#FSASF55 
XKTS SPAT eK XIS/NCHASY 
Secewe] Sor weeewccce 
AMFHSOXSTC8) XSTC3) 205 
MPEHECKSTCZIOXST CSI) 05 
AS1=CHPFHOXSTC2) SHAE 
AS22049 Femelle SHMFA 
| ASS=CAPFHEXSTC 1) 7H4E4 
} XXT5206 ' 
as 30 LSJMe le VOA3Y 
1 HAL =2CLe“AUCISUM) DW CAUCISUS)+A61) 
| 442 =CAGCISUY) ©A52) eC AUCISI4) 463) 
FSA=SQRTCHA1L/HA2) j 
FSOBSHYF AeA CISUY) O42 407 
30 KXISSKXTG¢FSA/FS8 
XXISSPATeOXXIG/NCHBY 
MUTT S MXIT LOKRKIZOMXISOKXIT OCA ZS eT ASO Ce XXIS HK XLS) SPAI 
TARITI=Z2 
TaRT T2230 
TWRITS=50 
TE (Ls er cel eANI ol Per selwRITI) ARITI CS e552 XXTLexXXT2eKXTS oxXxI 4elP 
TE CLs eEGe lL SeANDelI Fel del MRIT2) SRITICS 955) KX Lo XKT2 XXL 3eXX14 I? 
TF Clue Sel eAND eI Per aelWRIT3) GRITS (6955) XKTLoXXIT2eXXIT 3 qXX1 Gel? 
55 FIRGAT C1lBxXeeomoLlgl2el3el% I= FC4) AREmeewe QL EL Gee 2Xdo2Xe 
AeIPse_I4) 
RETIRN 
END 


ve 
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SUBROUTINE JFSIMSCANSS) 
DIMENSION S2SVCLOL eS2KERCLOL) oXST C9) 
COMMON YCOSeS3ETA2 
SIMMOW KITVC200) oXAITVC 200) eANSGE2S( 290) eSARC20269) 
SIMHI CAVKC1O0) oeSAVT C100) eSZTAS IETF AC oKCOCCeNDAV gL PUY Ge NS2 
COMMON AJC100) eo LSHAR Ae NCHSYe BSTANC100) eBBTAN2€100) ed ETAN2(100) 
SIMMIN FLAFANG DELTA JSAP GALAAL GS SATMA 
COMMON SBE TAeXXMe ICPIeSARCOD(513) 
CIMMON TDS oe XA eK eXT eT ANGE 2 eV Ce V%eISIGS oXL31SSe 31 GSeSMALSeDSS 
COMMON XSNC5) 6 CUE eERCeYVYV eX MoI TERAGSXSIOC8) oSxXSIIDICE) oV¥XS (8) 
TIMMON PSI Ze LPeSARCC513) eSAICOC513) 9 LPMs DE 
COMMON BEFANCSLIS PAZ TAMCSL3) oLUeh Ko XIT(0 299) oXUIC20020 9 XIX 
COMMON XROIND eG A2ZAA 03238 ¢C2CC8 
SIMYON AAMAS IBS ec CT re AS eB eCB8eDBeTSAUSC100) od GAIUS (ICO SNGAUS 
COMMON/SPCCESFACE 
P415301415392654 
THIS SUBROUTINE CALLED FROM IKFNEde 
JSE SIMPSINTS RULE « 
391 143 = 148 
XSTCIMO) = YXSCIMO) 
SJE = COSCIELTA) 
SJE = SINQJELTA) 
NS2 SAdDULD dav A FATTIR GF we 
VS 2207945. P M2 
VS21 = NS21 
NS24 = VS2-1 
S2GAP = (XSTC3)-xXSTC2) > /NS2 
JJ2 = CISCALFALSGAMYAD/COSCKST C52 ¢S444A) SESPATE 
D7CTSX/IAMHASHISCTSXCSOC/IAMMAG*IAFLACSCS = 2UJ 


apae 
fo 


o 


DP 2 182 = leNS22 
KS2 = KRSTC2*S25AP eC LS2-1) 
X40 = XS2eC3E 
K4AS = XS22KSTCA)#SJE 
XMASZ = XMASe#2 
ASO = XSTC4) #cDE 
AS$O02 = ASD#*2 
DsOX = DBGAP HXKO/SCOXHAS24ASI2) 2 7AT) 
T= C1S2edel) GI TO 3 
TF CTS2eEGeNS21) $0 TO 4 
TALL 32 (KS2e4NS529152) 
es G2 CALCULATES G2 WITHA XSI GIVENe 

£32 2 £xPCeansG2) 
Til ver de40) ANSG2SCIES2)=ANSS2 
S2KERCIS2) = ES2eowIK/UU2 
33 T9 2 

3 CONTINUE 
SA2KERCLD = DWIR/SARH CL o*xSTCS))3 
ANSG2SCTS2) =ALOGCSQRTC Le tKST(5))/UU2) 
33 T0 2 

i 4 SINTINUE 

S2KERQNS21) = Od OX/II2 
ANSS2S5¢01S2)=06 

2 CONTINUE 
§$28201) = De y 
09 10 JS2 = leNS2Ae2 

10 S2SR(J5202) = S$2S82(0552) 
LOOCS2KER GIS 2) Hee S2KTVUUS[SH1LI SAKE UCIS]+22 eS25AI7MS6 
TF CI veNE e090) GO TO 40 
SARC2C1) 200 
90 50 ISARC229NS292 
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ae eeeee 


an : te ET —_" eee ke. aut 
(tate ee SR emete T n MEME ORR RT Ra Oe ee nae mayen carrretbwerrelaiereer manne 


30 S2SICISARCISe5 HCS2STCISARC HL eSASUCL SARC #17) 
DD 36 ISARSSleNS21 
30 SART2CITSART) SS2SRCTS ARC) 
40 CINTINJE 
ANSS = S2SR¢NS21) 
RETIRN 
be Te) 
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ROUTINE IC2CSReSMeXCAGISIC) 


DIMENSION xSTC8) 

SIMMIN VCO ySIETA2 

COMMON KITMC 200) eX TT WC 200) eANSS2S( 200) oSARC2(200) 

SIMMIN CAV KCLOG) oCAV C100) oSETAIeTETAC aK CCCaNTAd gLPAMaNS2 
COMMON AUULIO) oe ISHAR ONCHSV SHS TANCLIO) eGBTANZC190) eFC TAN2ZC6100) 


CoM 
co 


MON FLAPANGITLT AS ISAP ALF ALSGANMA 
“ON SSETAK C44 TC2Le SARCIICS 13) 


COMMON IO UL a KA oXBoXleTANGerPoYCoVReISIGS oXL315Se81GSeSMALSeOSS 


CIMTON MSVG3D» SoD GE eV Vo KM oI FER AGSKSIOCS) oS XSLIICB) XS (8) 


COMMON PSI Zeal PeSARCC513) eSARCICSL De LPMeOE 

COMMON SETANC5132 e BE TAMCSL3) el uol Ko XTIC 200) oX dS (200) kd 
CIMVON KRIYVIGAZ44 93233 eC2C7 

SIMMCN BARA BABS eC TO Te AS OSI CSeIBeTSAUSCICO) pWSAISC1LOOI GNGAUS 
TIMMON/SPSAESP ACE 


30 1 TFN = 198 
LoxSTCIAND = YXS CIID 

XXL = XSTO4I*#SINGCJELTAD 

YYL = KSTC4dHCOSCIE_FA)D 

YY12 = YY¥te#2 
7 ISIC = 0 FIR WHINT 
> = 1 iv CAVITY 2 JFSiNS FIR FC5> AND IN CAVITY 
o 2 CALLED FROM FLINTL SOR Fld © 
5 2 FIV 12 3F Flare 

SR=Ce 

S¥=0e 

34SxXSTCLI #05 

244CS94°KSTC2) 

BaP LIs++1e 


BAMe SIH -XSTC3) 


Ae 


VSCXSTCBDOKST C32) 005 


HME SCXSTC32-XSTC3)) 0S 


341 


42OXSTCTI-LerdteS 


SPL4SOXSTO7) Led #05 


A351 
As2 


SUSMLERXSTCQ2) 4 C591) 
SOS4M1LAKSTO1LIIS C5714) 


ASS=H(SMimeXSTO3) I C9553 HD 


451 
As2 
as3 
312 
312 
313 
TF ¢ 
APL 
ar2 
a?s 
2aN 
090 


ne 
ao 


zas 


R35 
ia 
tt 
1? 
It 


242 F raKSTC 2) Sed 


ZCHPF Hele SAYER 
SC4PFaexSTC1) Ades 
sB44C434 


FBdP 1/3 


sS44F 734 


ISIC eNE 3) GO TO 20 
SOXCAPL OP MOKSTOLI KDA CKOAMKST C32) 
SXCaexSTC2) 


SSQrRTCAP1L/ AP 2) 


TINSE 

7 TSUMT1¢NCHBY 

CAJCTSI49 9311) #CAICTSI4) le? 
CAJCTSUM) 4312) ¢CAUCT SUM) ©8137 

CISL CotderS) RASCASCT SIM Obed eC ma UCTSUND©A51) 
CTSTCokGe m5) RIZCAUCTISUM) ©4527 CA CT SUM) A353) 
CISTI Colide=$) RA= AICISIMI ASII A CHASCISIM) th ed 
CEST Cer ae tS) RISCAICTS III SAS2IOCATCTISUMI ASS) 


SABSSaRTCRA/SRS? 


Sac 
xsi 
x74 
XPX 


Sea + SS aT 


TSteSQATCLoeeAUCIS IM e029 /545 
PS3H*AUCI SUM) O3n 

PsxSTPenx) 

P2EXPKIP 202 
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AVQ=EKPKP2CVY12 

VARIKIK FLA 2 

Rd l=vvyi/RV2 

TFcCISICeEdelL) RARZ1e/CXSIIPMKCA) 
TECISICeo£Ge2) RWRI10 
LECLSIDeZae3) RARSCLemSAC/APSI/( XSL F"KCA) 
IF CISICehae~5) RWRile 

IF CISICeE de eS) RaR=le 
SV=ESRHSAB*HAUIR 

SM=ESMOS ABSA We 

PAT=30141592654 

SR=SRePAT/NCHBY 
S4SS4eFAT/NINZY 

RETURN 

ENO 
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erereaergrereraed 


SIJBIDITINE FLINTLCYINT <CTRL) 

DIMENSION xSTCB> 

TIMMON YCSSeS3ETA2 

SSMMOW MIT EC200) ge KIT VO 2602 ¢8NS5286290) eSARC26209) 

COMMON CAVXC 100) CAVYC 100) e3SETAS OSETAC oXCCCoNCAV oLPMMyNSZ 
TIMHON ASVOLIO) pISAAR Ao NCAT oe BSTANCIIO) paBTAN2Z0100) pSETANZC100) 
COMMON FLAPANS DELTA JSAP eAL=AL GAMMA 

CIMMON SBETA oKX4 eI CAL ySARTIICS 13) 

SOMMON TOUL e XA eXBeXreTANGel Pa VCe VReISIGSeXL3L5 So 31 GSsSMALS OSS 
TOMMON XSNC8)¢ DAES ERCOVYV a XMe ITER AGSXSIOC 8) p@SXSIOOCE) oYXS 63) 
SCOMFON PSLZ eA SSARCCSL3D oSARCICSA3 I» LAME 

COMMON BETAN(523) oBETANC513) gL del PK oe XII C260) 9Xdd 6200) oe XOX 
TIMHIN KROIINDGA2K8O 93233 —9C20° 

COMMON AAAAS BABS oCCOLse ABs BSeCSe0Se TSAUSC(ICD) ed GAUSCLOOI®NGAUS 
SOMAON/SPOCVESP ASE 

SUBROUTINE FLIVTL CALCULATES THE INTEGRALS IN ©€1) 

LS4AR? = 9 FOR SHAR? woEeFIILS eo 

TS4AR2? FL FIV WIIVIEID wele =I LSe 

IF FOILS AAVE ROUNDED Lelee CATPFYSAEVGAUSS 

GIJAIRATIE 

QUADRATURE FORMULA TAN NOT 3£ USEDe SINCE SETA 

LS wl & SMIITH FINITLINe 

NOABY = NIMIER OF CHESYSHEV “GAUSS PIADVATIRZE CIVTROe POINTS. 
PAL = 30141532554 

Li CIC%leE2ed) GI TO 3 

30 70 1G2148 

KSTCI2) = XSNCT2) 

30 TO 12 

32 11 I4 = 198 

XSTCI4)2 = YXSCIHD 


CONTINUE 

INL = CXSTCL 10d #05 

ON2 = CKETCLIeLedeed 

All = CON2*xXSTC22)/9N1 
Al2 = CONQ-KXSTCSIIZING 
B55 = CXSTCLIOXSTC 2) 05d 
CABSSCKSTCQ22eKST CL IIH eS 
ASL = (BCS+1e2/CMBS 

A32 = CHBsr5exSTC39975435 
FCAS OXSTC32 -XSTC22)2 eS 


FOILS = CXSTCSP+XST C2 eS 

AGL = CFCLS*LeFCAS 

AG2 = CFCL59XSTC122 75245 
SPACES] = CKSTOBIACST C2997 .944 
READ LPMM FOR THE SZSUND ARCe 
LTECKCT RL eS Ze2) 50 TI 100 

TF (I veSEe2) GO TO 100 


SSPACE = CLeoKSTCLIISFLIAT CIP 4) 
FSPASTE = CSPATE/FLOATOCLPM=LOK) 
T2451 


KInIK = ole 

SPADESCSPACE 

32 20 TC43YSLeNCH3Y 

NCHINCHBY=ICHABY*) 

AJCICHSVIETISCKI[eN CHL eFaAT/C2eNC48Y)) 

MCSTSIVISASCI OAS) © 3N2 

TECLTIR Ast aet) 30 TI 483 

BECK CASK eo eKKSID GI TO 21 iy 
TE CTOMeGEeL PK) SPACE = FSPATE 

MIACK = XS4CK*SPACE 


At ant Abed en oN Pn et i epee a 


Lom = I194¢2 
33 T3 22 
KAKSI EXISTS 3Td KSICIQ4e1) AND XSICTIIWD 
21 CONTINUE 
toma = LO4e1 
BBTANCICHBY) = BETANCIOMD*CSETANCIOY)<“BETANCIOMAD) 
Ke Cx<SLexCal</SPASE 
S BSTAN IS JSZ0 FIR CHESYCAHEVHGAUSS INSTEAD OF SETANe 
32 TO 20 
433 X3TANCICM3Y) = SBETA 
2 BETAN FOR ITERAcEZel IS SPECIFIED IN OFSIM1e 
20 CIONTINJE 
120 CONTINUE 
TF (KIT e505) GI TI 5 
LF (KT T2®Leta05)} 39 TI 5 
LTECKCT RL ele 4) GO TI 4 
LF (KS Tinvelde3) 50 FI 3 
I? (KSTRLetSe2) G9 TO 2 
TS CLISAAR%*EGel)? 52 TI 10 
YINT = Oe 
33.119 ISJ4 = LeNC43Y 
B2C = CASCISUMPFALID AC AUCTSUS) OA12) 


oo 


LLG YINT = YINT *B3TANCISUM) e@SaRTCARS) 
YinT = yIvferal/vca3y 
39 TO 1900 


10 CINTIVSE 
* THIS IS THE CASE OF HANILING RVJED Lo ie 
VOF = 0 
XDA = De 
TALL JESIMICYINTeNO=eXCAd 
7 XDA LS DUMMY e BNLY JSZI FIR FCS) LNIKFENE We 
39 TO 1000 
S CINTINJE 
T$icse§ 
XTA=00 
TAL. LO2CSRoSY4oXCASISIC) 
YINT=SR 
3) TO 1003 
6 CONTINUE 
Tslc="s 
XCA=0e 
TAue LO2CSRoeS4oxXCAeI SIC) 
YINTsSR 
39 TO 1000 
2 CINTINJE 
KCA=0e 
2 KIA TS DUMMY ® 
TSic=2 
TALA [O2CSReS49XCAgI SIC) 
YINT=SS& 
$0 TO 1000 
3 CONTINJE 
eroerINTEGRAL FOR I3e 
AUOV) IS SALCJLATED aN) STORED 
YINT = Ge 
33 122 ISJ4 = leNCH3IY 
431 = LewasCI5u4) 
AB2 = CAUCISUMDFASL) #CASZ-ATCI SUM)? 
S3A32 = SQRTCAB2) 
ABC = AS1/SGAR2 


anes 
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prea ene eee 


See eee meteor 


120 


ever erararcreraer 


a 


38 


139 


135 


130 


152 


151 


230 
i?od 
131 


1930 


ve 


YINT = YINT#a3c 
YINT = YINT*PAI/NCHB3Yr 
30 TO 1000 


wmeeeINTIGUAL FIZ 14 


SINCE SETACN) BTWN TCT ANDO TEP ARE 
EXPECTED TI BE ALMAYS SMIIT4— SSE SAISS= 
CHESYSHEV QUADRATURE FOR4ULAe’ 

AJON) TS ALREADY CALSILATISOe 

Te TAIS IS THE FLI&ASS CASE =I2% BETAN2» 
JSE A CONSTANT FOR 3ETAN2e 

33TAN2 1S JSED FOR T4Evv-GAUSS INSTEAD OF BETAN2Z« 
CONTINUE 

TP CITERAeSEe2) $9 TI 150 

TF CIueSTe2) GI TO 141 
ZTA2 MIST 35 READ FIX THE FIRST Wve 
33°16) [C43Y = levoi3r 

S3TAW2 (1543Y) = SBETA2 

NS215NS2e1 

23 1¢5 TOCS1leNS21 

JETANZCIOC=SHETA2 

39 TO 131 

CIONTIVJE 

L?¢iJesie2) GI T3 131 

IOMM = 1 

XOHOK = XSTC2) 

OO 170 IC4BY = leNCASY 

X4ST = FCAS#AUCICmSy )+FC15 

TF OxC4CK eS e XSI) GI TO 151 

KOASK = XCACK * SPATE2 

TIM = To44e) 

30 TO 152 

CINTINVE 

TIMMA = ID44=1 

BSTANQCICA3Y) = SETANQCII44) 
LedBS TANS CIOMM) SSETAN2CITOMMAD DOCK KSI OX CHCKISSPACZT2 
TuMtIS Ady 

K4SL = FCAS@AUCIO4 = 4FC15 

MGRITECS 9250) ILMeBETANZCILM) eXKSI 
FIRMATCLSXe @Loegl 3g 2 Kot B3STANQS tg El 4e TeX geXKSLaeet loes7) 
CONTINUE 

SINTINJE 

YINT = Ge 

29 199 ISUM = leNCH3Y 

AZ1 = CSBTAV2CIS S49) +Fal dw ( Lesa CISI4)) 
AS2 = CAUCISUM)%AG1) #CAUCISIM) 4A42) 
SaA32 = SaRTCAB2) 

YINT = YIVT * A5S1/S94A32 

YINT = YINT*PAI/NCHar 

CONTINVJE 

RETIRN 

pe ie) 
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SUBROUTINE G2 (xS2eA52eTS2) 

DIMENSION MSTO63) eX12166290) oXI1225 0200 eX125SC290)exI24S0200) 
COMMON YCOCySBETA2 

SIMMOW XITHCQ206) eX ITNC 200) pANSG2S(0230) eSARC2(209) 

SIMMON CAVKC LOO eCAVFC LOO sSETASeZSETAC eXCCCoeNTAV gL PYM eNS2 
COMMON AJC1L00) pISHAU7 NCHS YO BSTANC 200) eSBTAN2€100) o8ETAN2(100) 
COMMON SLAPANG DELTAS IGAP gAW= Al BANA 

COMMON SBETAeXXMe ICPI ey SARCOI(513) 

COMMON [DJ eo XA XT eXT oT ANGeE eo VCoVReISIGSeKLIISSe SIGS eSMALS OSS 
COMMON XSNCS)9, CLEGERC eV VV okMo ITF ER AgSKSIOC 8) oSXSIQOCB) eo VXS (8) 
TOMMON PSI2eLPoSARCI513) eSARCOCS13S2 9 LPMe DE 

COMMON FETANCSIS SG ZETAMCSIS eI e uke KIT 200) eK UIC200) KIX 
COMMON XRIUND Gg AZAR 93238 eC2CO 

TIMMON AAAAGBISZeCCI To AB SAS eCS eI Be FSAIS(100) ed SAIS C1OOIeNGAUS 
COMMONS SPC/ESPACE 

DOMAINS TAIL UT AS 

SIMMONS CVTYL/CACLENS 3135S2 

TaIS SUBROJTINE ITS CALLED 3Y OFSIMS. 

T4LS SS3SRISTIVE CALIJLATES FUNCTION S20xS2) stig 

INCLUDES 121¢€KS2) TD T29€xS2.¢ 

KS2 [5S xXSI* AG2 IS FHE SILITIIV JF INTEGRALS © 

33 1 LaP=193 

KSTCTAPvErKSCIQS) 

Pal = 30141532654 

CCOCLSALOG CL a KSTCSII/C 2eePAl) 

SSG=SIWTCL ee XSTCS)? 
U2SCOSCALFALSEAMMA D/C TISCKST CS) GAMMA ZESPACE 

J22=UZe%2 

Y22L=1le/U22 

Ta OSSSS*CAV LENS J2¢ (Le SAVIEND 

TAUsOSALIS(T WO) 

Twl=SaRTCLembLe-XSTCG) -U22ZLPeCAVELEN) 

TAUMLTSALOSC TAD 

TE CUTAU CES 00) TAUWETAUEO 

TF Culadel ded) Thdd=fAddl 

IF (1JeSE047) GO TO 100 


ore L2L CXS oree, 


Tae SAVE CNTESRAATIIV AS T4ar IN 
SUBROUTINE CAVITY FIR GCCXSI) 


NOF = 3 
CALL OF SIMI CANS eNOF 9 XS2) 
xt21 = ANS 


TFT yeEGehO) XI2ZISCIS2) = XI22 


oer nl 22 0X8] ennne 


JSE THE SAMD SUZIIUTLVE L252 as 
JSEG IN CAVITY WITH I[SICH1e 


tsIc= 
CALL [C2CSRoSMoxXS2eISIC) 
xI22 5 Sk 


NOTE THAT SM IS OUM4Y VARIABLES 
TFT seb GehO) XI22SCIS2) = x122 


eeeaL23 (XSL) wees 


USE CHEBYCHEV=GAUSS QUADRATJRE FORMILA 
IN ZXACTLY SIMILAR WANNER T2 THAT IN 
SFSIM3 FIR I3e 

X123 = Oo 

SPCS = CXSTCLI OXSTC2)0 005 


cas OXSTC2=XSTCLI Ie 0§ 
A31 = (€BPT5 + 10945435 
A3S2 = (-8PCS © XSTC322/CMBS 


176 


ererenga 


oF 


a? 


a> 


« 


190 xI21i 


32 2 ISUM = LeNCA3Y 
Aal > LemasCISu4) 

aA2 = 

$442 2 SQtFtna2d 

F313 = MAL/SHA2 

F3AIS = CHAS @AUCISY4) e379 C5 eK S2 


2 XI23 = XT23¢FSIS/F3AL3 


XT23 = XLQ3ePAL/NCNSY 
I= (C1 seEGedG) XIZSSCIS2) = xXI23 


eee eal QG4ewone, 


JSE CHESYSAEVeGAUSS JJADRATIRE 
FORMULA 3Y ASSUMING THAT 
THE KEWVEL FINe C3 S4IIT4e 


HU = CXS2¢L0ed# CXS2=KSTC1)>) eC RSTL3)2 = KS2) 


avo = KS2-xSTC2) 

aw = SQRT(4Us4V? 

FACS = OXSTCSeXST C22) 6 05 
PMCS = OXSTO3)XSTC2)) 05 
AG1l = CFPCS+le/FMCS 

442 = CFPLS-KSTCLIISFICS 
XT24 = Qe 

33°10 ISU4 = le NOA3SY 
TOAL = AUCTSU4) AS] 

T2A2 = AVCISUM) #442 

SPP = SQRTCTPAL®TAA2) 


F4T = (CBSTANQCTSUM) FAT) © C1LeeAdC ISUM))/STP 
B3TaAW2 [S$ CHE3V=SAUSS VERSIIN SIR BETA IN T4E SECING ARCe 


FHA = FUCS#AdCTSUY) &=>C5"xKS2 
ST2 2 SQRTEL er AUCI SIA 202) 
FaB S FACS wST2e(3EFANACIS2 deal / 4a 


10 X12% = XT2¢¢CF TWF aad /F aA 


XI1241 = xL24#PAL/NC3Y 
SETAN2 IS USEO FOR SIMPSONTS RULEe 


KeG S ALOSCOXSTO3) =K529/0X52°xKST(2))) 
IS2 1S TRAVSFERRED T4R0UGT S2-ARGUMENTe 


X1292 = KLo*eCSZETANZ(T52) oP ALI 4a 
XT29 = xXI241¢x1242 
TFC ueEO049) KI24SCIS2) = xXI24 
32 TG 101 
= xT21S(1S2) 
xT22 = x1l22S8¢1S$2) 
xX123 = xI23S¢€TS2) 
XI29 = xI29S8¢IS2) 


LIL K3S2A > -KL2L/>alexnl22 
XSZBSCCCIMALOGCCOS CALF AL*+GAMMAD/COSCXSTIS) +GAMMAD/ESPACE) SPAL 
LAPSYSTUTSXA>AMIASSISCTSXCSOC/) AMMAGOLAF LACSIC(SOLA=1CCC 


xS2C = xS28xI23 
R523 2 -XT2A/P AT 
AS2 = CXS2A*XS2C64#KS29) ed 


Deen] 25 cee ccwwwewence 


S7L4ISCKSTC7) 41 od 205 
SMIHSCXST CTI 91 od ee 
ASIECSM THO KSTC2) 9/0 95291H) 
AS2=CSMLH@-XST(1)90/ 0 *5P1H) 
A5S32CSTINMKSTOC3I2/ Coa? 1H) 

XS2E2960 

a0 20 TSUM=1eNCHBY 

AGL F=CLoPAICISUN) DO CAUCISI4) +451) 
HA2 =CAICTSUMD OAS 2) eCAUCTSUMD A535) 
SHI2Z=zSARTCHAL/HA2) 


177 


CAUCTSUSD © ASLIFCAS2°AICI SIND) 


= 


328x 


= Se re nee enn ea te ere re ee 


NA3S@GPLHOAUCISUM) #5 41LHH=xXS2 
FSSSAL2/HA3 
20 KS2ESKS2EeF5 
XS2Z=PATEKSZE/NCHEY 
Teeel2§ cer cescecccces 
HMFHECKXSTC3D-XST(3)) 205 
HPPHECXSTC8) OXST(3)) #05 
AGLSC4HPF He XST(2) 2/04 
A522 ¢4P Fedele d/AVEH 
ASS C42 FAOMSTI sada 
XS2F x06 
33 30 TSJMS1eNCH3Y 
MAL SCLeeAUCISUM)@CAUCISU4)¢AS1) 
MAS SCAUETSI4>9 ©AG2) @CAUCIS IND ©A63) 
S412=2SGRTC4A1/ 482) 
SASS AVF eA ICT SUS) OFF 42X82 
Fszsnl2/443 
39 XS2FSXS2Fers 
MS 2FsIaTexS2F/NCA3ZY 
AG2SASG2 oH XS2EeKS2E) eTAUN end SPAT 
LF Ch set Ge 27 eANDolS2etde2) ARITE65 952) XIZLe XT 229X123 ex 24¢eIS2 
IF ClJyeTGe27eANDolS20fQ0elD) aRITE(C6052) KI2LeKL22eXT 23K 2491S 2 
TF Clue bGe27 eAND ei S205G050) WRITECG2 52) XI2leXL22eXl23eX1 240182 
52 FIRVATCLOX eg wemmTigl2ygL 3914 JE FCS) AW smwma ge Hl EL e7e2Xo2XKe 
A #IS220_14) 
RIZTURN 
END 


< 
<« 
prema Diaconate hee me nek erence At dalle he ih nineteen ce hanna iat ia tie isan eS oan ee 
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aaa 


SUSROUTINE CAVITY (xf YCC) 

T4IS SUBRIVTINE TS CALLE) FRIM IKFNEW FIR FCSde 

OIMENSION CKEXCLOOd sS<KEYCLODISANSTIC LOE) eSRI26100) sSICSIS C1009 
DIMENSION SIC81I4C1L O39 oe XSTCS) 

JIMENSTIV CAVKKC190) erAvYY(109) 


COMMON 
COMMON 
COMMON 
COMMON 
TI4M4CN 
SIMMON 
CI44ON 
COMMON 
TOMMON 
COMMON 
CIM4IN 
cTIMMON 


YCCCeSSETA2] 
MITMC200) eX ITNC200) eANSS2S6230) eSARC2¢200) 
CAVXO100) sCAVY( 100) oe SETAZ SETAC eXCCCeoNCAV pL P44 NS2 
AJC1L00) oIS4AR7 og VCHSV 933 TANC19O) eBSTANZC130) oZETANZ C100) 
FLAPANSDEL TAs JSAP SALT ALSGASMA 

SBETAgKKMe ICPIeSARCOI(513) 
LOJn XA KS eg Xl oTANGe De VCoVVoISIGS oKL31SS¢ 81 GS eSMALS OSS 
XSNCBD6 TRE ERC OV YY oXMe ITERAG SXSIOC8) oSXSITO0C8 2 o YXS (8) 
PST Ze LP oSARCU513) pSAWCICS 1352» LIM oDE 
BETANCS13)¢BETAMC513) oI dal >Ky XIT10 200) XII 200) XIX 
KRIIVO GA ZAK 93233 0C2255 
AAAASBBAB OCCT Se AS OBS CSeDBeTSAISE1002 od GAYS (100) 9NGAUS 


COMMONS SPOVESPACE 


TIMMONS 


TAJI/ UT AYU 


COMMONS CVTYL/CAVLEN, 313S2 
KZCC IS THE CAVITY END FDINT CALCILATED IN SUBe CAVITYs 


COEL = 
SJE. = 


COSCIELTA) 
SIVQDELTAD 


Pat = 30141592654 
39 1 WOA=198 
eSTCLIAd = YXSCL3AD 


scomM = 


SQRT C16 +KXST(53) 


COLFALIGCLe+KSTC EDI SC 2eeF Al) 
SSG=SQRT (Lee XST(6)) 


A2z5I5¢ 
3$222N20 
U22T S16 
Lw02S835 
TAUM OSA 
Tdi SS ar 
TAUMI=A 
I? (Ta 
1 (ura 
NOAV=80 
VOAVI=EN 
Cavs = 
Liave T 
SINS 
ao 2 Se 
Ata = 
RIAL PA 
T= (4M 
TF CKL4e 
Tricxsl 
TF (Ide 
NOF = 3 
caue OF 
ANS IS 
TT? (Ide 
69 TO 7 
ANS = A 
CONTINU 
To2]2cxst 


ALFALOSATIAIZTISOKST CS) +SAMHADS ESPACE 
*2 

422 

WCAVIENS JS2e¢L es CAVLIN) 

LOS( TWO) 

TOLeeCLeeKSTCS)-J220. eCAvLis) 
LOG¢TW1) 

Jdotded) TAUSZTAUWO 

Jetael) TAJd=Fasdl 


Cavell 

CXSTC2)-xKXSTCLII/NCAV 

HE CAST PUINS JF xSI = C SINDE THERE IS & 
ULARITY FOR SINGLE SPIRAL VORTEX MODELe 
“M = LeNCavi 

ST¢1) #Cavs* (4.491) 

RT OF JMEGA = SETA FAle 

eES0e1) GI TO 3 

EGeNCAV1) 350 TO 10 

Y SALCILATION, CALLING JFSiMte 

GEe3%) GO TO 75 


SIMICANS e VIE gC cA) 

A SOLUTION FIR ICACKCT)» KCI IS IDENTICAL TO xXCAe 
EQe27) ANSTICG4) = ANS 

6 

NSTiCSOM)D 

£ 

) CALCULATION 


LiCl ¥055039) $3 TI 77 


Tstc = 


1 


Oa a ee a en ee ee US 


| 


eetentnty mao mate eee 


TALS TO2CSReS4exCael SIC) 


> SVLY SR IS UTILIZED*= SM IS FOR RYIVTe 
TF C1 se EGe27) SII204K 4) = SR 
33. T2 73a 


77 S2 = SRI2QCKLM) 
T3 CONTINUE 
ComeenIC2 CxSI) CALCULATIINe= USE CHEBYSHEV-GAUSS 
: JJADRATURE FIRMA’ 
3PCS = CXSTOC12+xXSTC2)) #05 
CMBS = OxSTC2-KSTCLII 005 


A3S1 = COPI5¢1e9/C0485 
A3S2 = CHB2CSeXSTCSII/OMBS 
241 = xCawxST2) 
EK2 = OXCAtLed eC XCAmxSTC1)) @CxXCAKSTC3)) 
243 = SGRTCIKI/ERK2) 
f*33 2 C435 ef <3 
TF ClyeGEe34) 69 TO 30 
SIf3 = De 
39 3 FSUM = LeNCNSY 
TILSCASCIS IMD OAS1) CAS 2-AICI SIA) 
SEul = SGRTICEY) 
TF3 5 CLeeAICISUMI ISI U1 
EF3a CMIS*AUCI SUM) O3°CS<KTA 
§ SIc 


SIC3ePal/NCn3y 
SIC3 = SISS*ALOGECXSTC2-KXCAISEXCAKSTC1))2 EK 
L2 Clue D027) SICSIZCCL4) = SIT3 
$0 TO 381 
30 SIC3 = SISSISCKLM) 
31 CONTINUE 
ermmn [COCXS1) ower, 
JSE CAHEZSYSAEVOGAUSS AJADRATIZE FIWILA 
IN THE SAME 4ANNIR AS THAT FOR I4 IN 
QFSIY36 
L&CIveGEet4) SO TI 32 
FICH = CKXSTCZI ©KSTC29) #05 
FMc5 = OXSTCSH-XST C2) 0% 05 
AGi = CFPCS*¢Led/FMCS 
492 = CF 2T5=xKSTC1)7F4C5 
SIC% = Ge 
332 7 ISJM= 1LeNCHBY 
RA = CZBTANSCISUM) > ATI #CLeeAICI SU)? 


3 = SISS*CEF SHEE 3ZeSART CL om AICISI4) e820 /E= 358 


eaqgaeaed 


33 CASCISIM) ©AS1) © CAICISI4I)D ©0392) 
329 = SaRF CRB) 
RE = AA/SR3B 
‘ WD Ss FMCSeALICISIM) +7 255 =KCA 
7 SIC&® = SIC4*#RE/RO 


3I1CO = SISSePAI/NEAZY 
Til sefGed7?) SIC4I4C<KL4I= SICS 
32 T3 33 
32 SICH = SIT41 44.4) 
33 CONTINUE 
LF C1 se EG e 27 eAND KL Gel del) AIITI C5955) ANSoSMeSICSeSIC4eKLM 
' TF (Tse EGe27eANDeKL4e7GeS0) eRITEC 6955) ANSeoS%eSICSeSIC4 eo KL 
BF CL setae 27 eANoMuitelr Ged G)? ARITE( S953) ANSeSAeSTCSeSICH oe KLM 
53 FIRYMAT C1LIXe to mHeTlel2el3el% OF CAVITY AREmm= He GCE lL4e 7a 2k o2Xe 
Ae<L mse, 14) 
> 4SOKSID = L4EK3) ALREADY CALCULATED e 
WUU2 = COSCALFALOGAMYAD/COSEXST CS 2 ¢SA4MMAD/SESPACE 


errr eerr errr rrr 


180 


SPLHSCXSTC7 el eod 05 
S4LISCKSTC7) 92 0) 05 
ASLICS4YLHKSTO2)940°67 1H) 
AS2=26SMINUOMSTCLI DS CH SPL HD 
AS3S0S4T LAO KSTOC399/ C9521 4) 
SICS=06 
03 30 ISJMS1leNCH3Y 
Ra So ClLetAUCISuY)) @CAUCISUM) ©A51) 
33 SCAVCISU4) A523 eC AUCIS IM) ©A53) 
SRASSSARTCRAS3B) 
ROS=GPLH®AJUCISUM) OGMLAMXCA 
FZSSRAI/RS 

30 SIC5SSSICS+=5 
SICSSPAL*SICS/NCHOY 

were ll bewewoeorce 
AVFHSOXSTOC32 oKSTC3)) 005 
42FAsOnSTO3Z2 KSTC320 005 
ASITSCAPFHOXSTC2)) /H4MER 
AZ2=04>F te] osAMEA 
AGSSCHPFMAXSTCLIIS HEH 
S1C6=0e 
30 40 ITSUM=1eNCH3Y 
Ra SCLerAUCI SUM) @(CAUCI SUS) OAS) 
a3 SCAUCIS S42 6AS52) eCASCISI92 4863) 
SRABZSQRTCRA/RB) 
ROSHMA 4 eATCISIM) O4P 74" XTA 
FESSRAB/RC 

30 SICS=SICE+FS 
SICS=PAI*SISO/NCHSY 
SS ot CMANSSPATHSROCSSOCLSALOSCUU2)/ PAI ©SICS 
Lesloasralry/i<3 


SCHGCeTAUdeCHSICS*+SIT6) /PAL/EKS 
32 TO 25 
3 GO = ZETAB +> AT 
32 TO 25 
10 SCSBETAC+PAT 
2 BLTAX ANG 3ETACC 30DY ANSLES AT 3 AND C) MUST SPEICIEED IN COMYONe 
25 CONTINJSE 


XxXxS = XCA*COEL 
YYT = XCA@KMSTC4I#SOTL 
YYT2 = YYTee2 
Ku = STOOD eCDEL 
KQJ2 5 XK Jee? 
XYB = VYT2*KxKU2 
BWsDK = DEA MAXKS/CKYZ@PAT) 
cece = cos¢6c) 
$30 = SINCS2) 
C#C = JeOK/SCGM 
SKEXCKL 4D = CSCHCFC 
SKEY CK 4) = SSCHCFS 
CONTINUE 
CAVXX(1I=06 
CAVYYC1) 206 
39 15 ICAV=39NCAV1 92 
TAYVUXCI CAV) = CAVXXCIZTAVO2) MCAVSHUCKEXCICAV=2) Hhee 
LOKEKCICAVMLIOCKEXCISAV IIS Se 
15 CAaVYYCICAV) = CAavry([cave2) 
LOCAVS#CSKZYCICAVH2Z eH ee SKEVCISAVR LD OSKEYCICAV) B/30 
Tel uef 2027) 30 TI 109 
60 TO 101 


130 


1j2 


191 


39 102 ICAVSleNCAVIe2 
CAV CIT SAVISTAVXXCT CA YD 
SAVY CICAVISTAVYYCICAY) 
KIOCOSSAUXCNTAVI) 
YCCCZCAVYCNCAVI) 
CONTINJSE 
XCCHCAVXXENCAVID 
YOCSCAVYYCNCAVI)D 
RZTIRN 

ENO 


182 


= = PU atalg sea CS cn iam iid. PaaS ort pine Segui 


SUBROITINE AMINT (SieSMe 412? 
JEMENSION CSTCR) 
CIMYON YCCCsSBETA2 
CIMMON KITWCZ[9G) XI FNC 200) eANSS2SC€ 230) eSARC2(200) 
COMMON CAVX(100) oCAVYC100) eRETASeSETAC oXCCCeNCAV eLPMMeNS2 
TIMMON AUCLOO sISHAR2 oe NCHBYe SBTANC1I0) pSETAN2ZC100) eSETAN2 (100) 
SIMIAN FLAP ANSIELTASISAP OAL = AL SSAYMMA 
TOM4ON SBETAGXX4eI COT eSARC0I(513) 
COMMON LOS eXAeXSe XT eFANGe Ero Ce V%eIZSIGS g@KLITSSe 3LGS eSMALS DSS 
COMMON XSNCB)¢ CLE GERCOYYY oe XMe ITERAe SXSIOC8) eSXSIOOCB) es YXS (8) 
CIOMMIN PSIZoL> oSAVWCC3513) eSAVWTIC5 15) eo L OMS DE 
COMMON BETAN(513)¢ BZ TAMC513) oIUelLP Keo XITIC200) eX JIC 200) 9XIX 
SIMMON XRIINDGA2AA 03 233eC2Cz 
ZTIMMON AAMAGB333 95555989053 09530939 TSAUSC100) oe GAUS(1L00) eNGAUS 
SOMMONSSPO/ESPACE 
Pal 5 30141592654 
IF (ISP elae0)? SO TO 10 
32 12 15 = 198 

12 xST¢ls)? = xSncts) 


33 TO 12 
19 399 1 IS = 148 
1 xSTCIS) = yxScIs) 
Li CINTINJE 


XX1 = XSTCS)@SINCDELTA) 
YY1 = XSTC4) #COSCIETA) 
YY12 = YY¥1le82 


235 = CXSTC2)=XSTC1) 805 
BTS = GKSTULIOXSTOC2))2 #05 
A321 = (BC5¢1.6)/C85 

A32 = Cest5exXSTC300/535 


BM15 = CXSTEL led eed 

2915 = CXSTCLI Led 205 

ALL = (€9415=xXSTC2) 273215 

Al2 = CAML5"XSTC3) 9735715 

FICS = CXSTCS*KST C2) #65 

EMCS = CXSTC3I<XSTC2) 205 

S>LI=CKSTOC7) 1 0) eS 

SMIH=ECKSTO7 Led #05 

MMF 4SCXSTCS)9KST C3) #05 

APEASCXSTC 32 *XSTO3)) eS 

ASLICSM LMA KSTC2)0 976-5914) 

AS22CS4 1G KSTILI IS (e529 18) 

A335°03S41R-KSTC3) 04025914) 

ASLSCHFF rte KSTC2) 0/44 F4 

ASQ CHP FHel eS HMFS 

AS3S= CAPE HOKSTCIIISN ATS 

AGL S&S CEPC5+1L ed /FACS 

A42 = CFPCS=xXSTC1I/F4C5 

T= (4laelae5s) 33 TID 5 

IF (MIQeESe5) GO TO 5 

T2C4120EG0e4) 69 TI 4 

TF (MIQeEQe3) GO TO 3 

TF (M0 aeEGe2) GI TO 2 

AJCI) ARE ALREADY CALCULATED IN SI3STIITINE 
TEINTLE ANDO STIRRED IN COMMIN AREA’ 

S%=06 

S4=06 

33 20 TSU“ = LeNCHBY 

GR1 = LesAUCTSU4) § 

3Y1 5 CAUCTSUMIOASIIOCAS2Z©AICIS IY)? 


ener 


183 


a emenemncenarmeemeemnmnemmemesncccmnasmesncsceimee es z x 


Ss¥1 = SGaireor2) 


FES = GxXI/SGY1 
FeL = CB8Saas LSM) 6325 
Fx2 = FxX1-XxA1 


FX22=5x2Qae2 
Fx3S = Fx22¢YY12 
FF3L = Fx2¢7x3 
TFS2 2 VYLCFXS 
SR = SReFFT3eFF 31 
= SOFT SeFF S32 
SR = SR*ePAI/NCHAY 
S4 = SM@PAI/NCHBY 
32 TC 1000 
SINTINGE 
S206 
S4=206 
32 50 SIMS LeNCASY 
GX1 SoC LeeAUCISUM)) eCAUCISUM) *AS51) 
SK2  SCAICISU490A52) CAS CISI4) AES) 
FESSSIRT(SKI/GK2) 
XSPF-SPLHeA CI SUM e541 dt 
KILL ECS > ony 
XPXL22xXFx1lee2 
JITI=ANXOKLZ¢VYI2 
CESLYSXP XLS BOT2 
FEs2zvv1s3dT2 
SRISKe FET ReFS 3] 
$0 SHSSMeT Ee 3er= $9 
StSSVeral/vcnsay 
SMSSMHOATSNCHBY 
32 T2 1000 
6 CONTINUE 
Si=960 
S4=be 
33 °7C TSIMZ1eNCH3Y 
BKL SOL e OS CISIMD DH CAUCTS IMO ASL) 
GX2  =CASCISUM) AG 2) eCAUCIS IMI *AE2) 
FS ZsSaRTcsKls35x2?) 
XSPSdI MAI CISUY) ©4974 
MIXLEMSPea Kk) 
K°XL2=K9X]Lee2 
JOTZ=KXPX12°¢YY12 
FE3LEK72K1/39T2 
FEs2zYY1/80T2 
SRsSROO CHA AI} 
70 S¥SSMeFF BFF 30 
q SRSSKRePAT/NGHBY 
S¥sSMe2alsyucnay 
33 TO 1000 
2 SINTINJE 
TF CISHARPSEQel)? GO TO 100 
TSAHAR> = 1 MEANS Treat T4E OTL 4AS 2IUNDED Lele 
S32 THAT THE SIMPSONTS RULE IS USEDe 
LSMAR? = 0 4EANS THAT THE FIL HAS S4AR? Loin 
S32 THAT CHE3YSHEV GASSS FII4LA CAN 3E JSED AS Z3E-Ode 
se =0 
s4 = 0 
go 39 ISuM = LeNCHBY 
STL = avcrSuMpeary 
STi2 = aAuUCLSUMPFAal2 


um 


arereanes 


FKL = SBTANCISUM) #©SARTCST1II/ST12) 
ANL = SPLSMAUCIS UT ORAL SK KL 
JNI2 = UNLee? 
SNES = UNI2¢YY12 
FK11 = UN2I/UN13 
F412 = YYL/UNLS 
SR = SR*EKLHF M11 
350 SM = SM*FKISFKI2 
Si = SXRePAL/NCABY 
SM = SmMePAI/NCHSY 
33 TO 1000 


130 CONTINUE 
2 THIS IS T45 CASE THAT THE F3IL 4AS VOUNDED Lele 
NOF = 2 
XTA = De 
TALL JFSIVICSIeNIF ok sa) 
2 KSA LS OUMMY=--2-O0NLY USED FOR 7(5) IN OXFNEWe 
NIF=2 
CALL JESTALCSYe NDF oXTA) 
33 TO 10C0 
3 TINTINSE 
JSE CHE BYSHEV=GAUSS FIRMULA SINCE BITA 
IN TALS RESION IS S43ITHe 
33TAw2 CISUM) ARE ALREADY CALCULATED AT TFINTLTe 


Grane 


3% = Oe 

S¥ = De 

33°30 ISUM = LyNCHBYr 

FZSL = CISTAVACISI4) FAL) we CLeehulI sud)? 
PSM = CAUCISUM) ©AG1) #CATCISIM) A842) 
SaPsm s S2atcrs4) 

FRG = PSL/SQPSM 

PSN = FMCSHASCIS IMD FF 2T5MmKK1 

PSN2 PSVe#2 


F54] 2 PSN/CPSN2¢YY12) 
F542 = VYLACPSN2¢YY12) 
SR = SRFFGRFF 4] 
S4 =S¥erF awk Fae 
CONTINUE 
SR = SR ePaAlL/NCAZY 
S4¥ = S4*PAI/NCHBY 
30 TC 1060 
® SINTINSE 
. XCA IS DUMMY, ONLY JSED FOR I52 IN 7 CE) 
KOA = Oo 
Tsic = G 
z SISIIITINE 102 TS ALSDd JSER IN FlSde 
SALL TOSS oS4oXCAei5ic) 
1936 RETIRN 
mvo 


us 
o 
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AD@A081 832 TETRA TECH INC PASADENA CA F/6 9/2 
COMPUTER PROGRAMS FoR CALCULATING PARTIALLY CAVITATING BLUNT TR-=-ETC(U) 
JAN BO S MAEKAWAr © FURUYA Neon s-Tere4a3e 
UNCLASSIFIED TETRAT=TC-3284-02 


ae 
onae 


eee ial haa ti ih cart een ne RMR ee Te een eee ne eter nae re leee 


SUBROUTINE SHAPECK of e3ETASIS112) 

COMMON/ FREE SAV/XFREE 29 YFREES 

COMMONS IPPER/AZAAT ST Z2IBJ oc2=CIoAAAA pSBBBUSTC2CJ eA Be SOU e CBU e2 BU 

COMMON YCDC»SBETA2 \ 

COMMON XITN( 200) oXITN( 200) 9ANSG2SC230) oSARC2 (200) : 

CIM4ON CAVXC100) oSAVYC100293ETASeSETAC oXCCCONCAVOLPYMONSZ 

COMMON AUCLOG) o ISAAR 29 VCHBY oB3TANGII0) oBBTANZC 100) o3ZTANZ C100) 

COMMON FLAPANS DELT As 3A> eAL TALS SASHA 

COMMON SBETAeXKMgICPI 9 SARCOD(S13) 

CIMMON LOL eXA oXSeXloTANG = 2 oVCoVReJFIGSeXL9135$0 31 GS oSMALSoDSS . 

COMMON XSNCB)» CLEERCo VVY oXMe ITERAgSXSIOCB) oSXSIOOCB) o XS (8) 

COMMON SIZ oL2eSARCCS1 3) oSAICIC513) 9 LPM 9DE 

COMMON BETAN(513) gBETAMCSL3D of do LPKy XIIC200) oXJJC200) oXOK 

COMMON XRIINDSAZAA932330C205 

COMMON AAMMe 63830 CCone ABs 339CS 028s TSAISC 100) 9d SAIS (100) 9 NGAUS 
SOMMON/SPC/ESPACE 
PAI 230141592583 

X2SKe02 : 

X3axee3 q 

XSZSQRT(X) 

XdsxexsS 

KE REEQSKEREEC HH? 

XE REESE XFREECOS3 

KF REESSSGRT(XFREES) 

XF REEH® XFREEC#*XFREES 

X22202e02 

X23502"#3 

X2SESaRTCe2) 

K2HE e2*X2S ; 

XB23 08882 

X83=08003 

XB S=SQRT C08) 

XBH=XSS 203 


WE MUST CHECK TD SEE IF d= ARE 52INS TO CALCUATE THE TOP PART i 
2&2 THE ZITTOM PARTe IF FI> HE TRANSFER FI 2NO AAL= JF RIJTIVE® i 
TS112 = 3 IS JSED FOR CALCULATIONS OF UPPER FOIL PROFILE 


qpacgd 


TFCISII20e&90e1) SO TI 30 
LF CIS1T2e5 203) G2 TI 30 


IF (XeLEoee2) 60 TO 15 
LF OxXebEee3) SO TI 29 
IF OXeGTee8) GO TO 25 


15 YOAQAAe XO3238exX2C 25 5exK3 
YIXZAZARS329BH 20 OX OCACCeSerk2 
SETASATANCYIX) 
GO TO 60 

20 VEARAA® (Hol3etXBo/So#XHoseeX2) ©8983 exeCCCCexs 
YOXSAAAAt( He SetBel3etLeSexSoSoeX) *3BBBeS#CCOC/KS 
SETASATANCYOX) 
39 TO 60 

25 V2ARC3IZ eXOCReKIeNHexs 
VIKZHSOVS oe TS eK See 3 eK2 i 
BE TAZATANCYOX) if 
62 TO $0 


T4135 2ND 4ALF JF THE RIIJTINE IS =O% SALCULATING THE JPFE8 HALF 


186 


pan tears 


Ma el aay ote 


30 If CIS1I2eEGe5) GO TO 70 
LT? Oxf UEESeEee2)? 53 TD 33 
IF (XFREES eLEee8) GI TO $6 
IF OxF REET eG Tee8) GI TI $5 
70 CONTINUE 
LF (xeucEee2)? 50 TI 35 
TF €xXeLEee3)? 60 TD 30 
TS CXeGTee8) 60 TO $5 


35 IF CISI12eE2e3) GO TI BO 
RISVFREECAZAAUOXFUTES 93 V™HBBIOKFVEE QV CACCUEKEREES 
32 TO 81 

30 1206 

31 CONTIVIE 
I® (xeGTee2) GO TO 48 
VOAZAAUOXOZS2BIVEKISI ASTI OKI eR 
YOXSAZAAUS[e MH VIBU OK OTe#C2ClUeX2 
SETASATANCYOX) =FA] 
32 TO $0 

90 Y2SA2RAUte2*S2B83UHX22¢C2CCU*eX23eR1 
VASVS-MAMASH CH eS Sate 2*SelSotKItwe dork 22) -B3H3U8 o2 -CClOCUeK2S 
TE CISLI2ZcEGe3) R2=00 
Le OX esT eed) GI TO 45 
VSAAAAUO (40/30 2X%B of See XHo ho tX2) 3983 9UHKHCOCCIeKSeR2 
YINZAAAAUO (Rel 3etBe/SotleS#KS-Betk OFFSBIUGOS PSOIOI/KS 
JE VASA TANCYOX) -PAal 
32 TO 66 

95 V3SABIJO3ES* eB oCIUHK3 2498 UeK33e22 
R3SV3“ABUMKI BUC eB oC BUEXB2-08IeXE3 
IF CISLL2eZ Qed) R3=Je 

86 YSABUCK BURKX+CSUSK2+03UeXSeRS 
VIN=ZBSF eZee lS BlexesoeIBJIex2 
SETA SATANCYOX)-Pal 
30 TO $6 

30 IF CISAL{2ceLAe3) SI TI 30 
R2=YFREECAAAAU# CH o/ Se @KXEREEC Be / Sot XE REEH HS oP X= REES) MBIBEU ORE REEC 
1 *CCOCSJex=RZEs 
69 TO 91 

30 22=06 

9% CONTINUE 
IF (xeSTee3) SO TO 35 
VSAAAAS*# (Ho SetX OR ol SoeXH~hoeK2) 3339 eX+CCOCIeKSoR2 
YOMTAMAAUREBeS 30OG al Set lLaS*KXSobe eX dO SBEBUS eS eCCOCUSKS 
ZETASATANCYOX) PAL 
G2 T) 60 

3§ IF CISI1L2e£20e3) GI FI 100 
RESYFAEECHAZIOBS JOR VEEC MCHIORFUET Se IBIMKEL EES 
39 TO 101 

130 %35=0e 

101 CONTINUE 
32 TO 46 


50 RETIRN 
ENO 


ve 
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Se Spat ait SPE, Nee SBS NO es phe meses 


tine 


SJBROUTINE XCYCCXCBe VOB eCXe SY) 
CIMMONSUPFERSAZAAU SS 2ZAUV eC2CC Ie AAAAS eABABUSCCSC J eABVe SHU CEUGD BU 
xX<3CX : 
XK2SxKe0e2 : 
AC3SKC ees 
XKS=SQRTCXK) 
XCHENCOXKS 
TPz0 
IF (CxXeLEee2) GO TO 3 
IF (CxXelLEee3d) 63 TO § 
TE (CxXeGTeced) GO TO § 
3 FISAZAAU OKC OB 2BBUOK< 2eC2Co Sent 
F2TAZAAUSZ o@ B2BBUS XKeSeeC2CCUeKxK2 ; 
FSxsxKeCx 4 
TXKZFLOCFESSE 2°CY) 
O1sF2 
J2zACIL OF Sel Vee BIRR see CIQSo Jenn) /OLee?2 
25 xK29102 : 
AIVSE UK SOF KC ‘3 
XKZXK-OlV 
T8siPol 
ZTZAISCOIVIXC) 
TCC ZeLEeeS OOOO eA el IP else29)) $0 TO 6 
32 T9 3 
& FISAAAAUECH eS 3 et UK OZ e/ SeeXKtoHeeXK2) +BBBAUEKKSCOCCCUXKKS 
FOSAAMAVE Che 3S otBelSotleSeX<SeBorXK) *ZABBUSSOCSO CI ee S/XKS 
FS2xK-CX 
FuCSFLOCF3SEQ<CY) 
322F2 
D2ZCOLMFSCCAAAAUEE Bel SotLeSteS/XKS930) eCOOCU ee Se oS /KKHDI/DI1 802 
3*x<391¢02 ' 
SI VSFXK SOF XK 
1>s1 Pe) 
ZEABSCDIVAKK) 
T2 CC ZenrEeeDOCOOLIORO KT %EGe20)) SD TOS 
33 TO 4 
5S SL SABIORBU eK Ke CBUEK<2eDBUEX<S 
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WAKE CLOSED AT INFINITY 


Figure 2 
Double Wake Model for PCASLDW 
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